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interpretation  of  the  SPIRE  data  and  the  interpretation  and  analysis  of  high 
resolution  spectral  earthlimb  emission  data  that  are  expected  to  be  obtained 
in  forthcoming  DNA  sponsored  high  resolution  spectral  earthlimb  emission 
experiments. 

Analysis  and  interpretation  of  spectral  earthlimb  emission  data  obtained  in 
'the  SPIRE  experiment  was  aimed  at  quantifying  the  role  of  C0£  in  the 
formation  of  the  non-LTE  earthlimb  at  wavelengths  15,  4.3  and  2.7  ym.  Non-LTE 
mechanisms  for  CO?  emission  in  these  wavelength  regions  were  verified  and  in 
some  cases  upgraded  by  comparison  between  predictions  and  the  SPIRE  data. 
Implications  of  this  analysis  for  the  gravity  wave-induced  horizontal 
structure  in  the  15-ym  earthlimb  and  for  improving  the  current  High  Altitude 
Radiance  Model  were  also  presented. 

Analysis  of  the  15-ym  SPIRE  data  confirmed  that  two  mechanisms  were  required 
to  be  incorporated  into  the  pre-SPIRE  model  to  improve  it,  namely  (1)  the 
detailed  altitude  dependence  of  the  line  and  band  shape  for  the  15-ym 
radiative  transfer  calculation,  and  (2)  the  rate  constant  for  quenching  the 
CO2  (010)  by  atomic  oxygen  was  reguired  to  be  considerably  larger  at 
mesopheric  temperatures  than  T-l/3  extrapolation  from  high  temperature 
laboratory  data  would  predict.  Three  mechanisms  were  found  to  dominate  the 
sunlit  4.3-ym  earthlimb;  these  are  (1)  absorption  and  subsequent  reemission  by 
CO2,  primarily  in  the  001-000  fundamental  band  of  the  major  626  CO2 
isotopic  specie  at  altitudes  above  105  km  and  primarily  by  'hot  bands'  and  the 
bands  of  minor  isotopic  species  at  lower  altitudes,  (2)  the  fluorescence  of 
sunlight  at  2.7  ym,  and  to  a  lesser  extent,  2.0  and  1.6  ym  sunlight  in  the 
altitude  range  below  105  km,  and  (3)  the  transfer  of  excitation  from  oxygen 
singlet  D  to  N?  vibration  via  a  quenching  collision,  0('D)  +  N2  is 
important  in  trie  55  to  65  km  altitude  range.  Analysis  of  the  night  time 
4.3-ym  data  indicated  that  (1)  the  transfer  of  excitation  from  0H(V  >  0)  to 
N2(V»11  is  required  to  model  the  data  in  the  75  to  90  km  altitude  region, 
and  (2)  a  second  mechanism  is  required  to  model  the  data  at  higher  altitudes. 
Direct  absorption  of  sunlight  and  a  "Chapman-like"  mechanism  were  found  to  be 
Inadequate  for  this  purpose.  Our  initial  analysis  of  the  sunlit  2.7-ym  SPIRE 
data  stimulated  a  coordinated  effort  within  the  AFGL/DNA  community  to 
understand  this  aspect  of  the  SPIRE  data.  This  culminated  in  a  combined  and 
complimentary  effort  by  R.  Sharma,  T.  Degges  and  ourselves  in  arriving  at  an 
understanding  of  the  sunlit  2.7-ym  SPIRE  data.  The  major  mechanisms  were 
determined  by  this  triad  to  be  (1)  the  fluorescence  of  2.0  and  1.6-ym  sunlight 
by  C02,  (2)  the  resonance  scattering  of  2.7-ym  sunlight  by  H20,  and  (3) 
absorption  and  reemission  of  2.7-ym  sunlight  by  CO2. 
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20.  ABSTRACT  (Continued) 

The  laboratory  effort  involved  studies  of  the  fluorescence  of  CO?  at  4.3  pm 
resulting  from  initial  excitation  of  the  101  band  due  to  absorption  of  laser 
light  near  2.7  pm.  By  utilizing  a  pulsed  laser  source  at  2.7  pm  we  have  been 
able  to  observe  the  resulting  fluorescence  at  4.3  pm  and  thereby  provide 
direct  laboratory  confirmation  of  this  mechanism,  which  is  a  major  contributor 
to  the  earthlimb  radiance  from  CO2  near  80  km  altitude.  Several  experiments 
designed  to  verify  the  branching  ratio  for  radiation  at  4.3  pm  versus 
reemission  at  2.7  pm  are  described  in  this  report.  Experimental  difficulties 
pertinent  to  this  measurement  are  discussed.  We  also  describe  a  measurement 
of  the  relative  rate  constant  for  collisional  transfer  of  V3  vibrational 
excitation  from  the  species  12C02  and  *3C02  to  N?  that  is  derived 
from  measurements  of  the  rate  of  decay  of  the  4.3-pm  emissions  for  13cc>2 
and  12co2  in  the  presence  of  N2.  Additional  possibilities  for  enhanced 
understanding  of  the  mechanisms  related  to  CO2  emission  which  would  involve 
measurements  of  several  LWIR  bands  of  CO2  are  discussed. 
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Section  1 

INTRODUCTION  AND  SUMMARY 

1.1  GENERAL  INTRODUCTION 

In  this  report  we  document  progress  we  have  made  in  two  major  subject  areas. 
The  first  of  these  includes  the  analysis  and  interpretation  of  spectral 
earthlimb  emission  data  obtained  in  the  DNA  sponsored  SPIRE  experiment,  and 
the  implications  of  this  analysis  for  application  to  cases  of  interest  to  the 
AFGL/DNA  community.  The  second  subject  area  involves  a  supporting  laboratory 
experiment  to  obtain  data  pertinent  to  understanding  the  nuclear  effect  of 
fluorescence  of  fireball  light.  The  laboratory  experiment  also  directly 
supports  the  analysis  and  interpretation  of  the  SPIRE  data  and  the 
interpretation  and  analysis  of  high  resolution  spectral  earthlimb  emission 
data  that  are  expected  to  be  obtained  in  forthcoming  DNA  sponsored 
experiments,  such  as  the  SPIRIT  and  related  experiments. 

1.2  SUMMARY  OF  THE  SPIRE  DATA  EVALUATION  EFFORT 

Analysis  and  interpretation  of  spectral  earthlimb  emission  data  obtained  in 
the  SPIRE  experiment  was  the  principal  thrust  of  this  effort.  The  analysis 
was  aimed  at  quantifying  the  role  of  C02  in  the  formation  of  the  non-LTE 
earthlimb  at  wavelengths  15,  4.3  and  2.7  um.  Non-LTE  mechanisms  for  C02 
emission  in  these  wavelength  regions  were  verified  and  in  some  cases  upgraded, 
by  comparison  between  predictions  and  the  SPIRE  data.  Implications  of  this 
analysis  for  the  gravity  wave-induced  horizontal  structure  in  the  15  um 
earthlimb  were  also  presented. 

Analysis  of  the  15  ym  SPIRE  data  confirmed  that  two  mechanisms  were  required 
to  be  incorporated  into  the  pre-SPIRE  model  to  improve  it,  namely  (1)  the 
detailed  altitude  dependence  of  the  line  and  band  shape  for  the  15  um 
radiative  transfer  calculation,  and  (2)  the  rate  constant  for  quenching  the 
C02  (010)  by  atomic  oxygen  was  required  to  be  considerably  larger  at 
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mesospheric  temperatures  than  T-^  extrapolation  from  high  temperature 
laboratory  data  would  predict.  Three  mechanism  were  found  to  dominate  the 
sunlit  4.3  urn  earthlimb;  these  are  (1)  absorption  and  subsequent  reemission  by 
C02,  primarily  in  the  001-000  fundamental  band  of  the  major  626  CO2 
isotopic  specie  at  altitudes  above  105  km  and  primarily  by  weak  'hot  bands' 
and  minor  isotopic  species  at  lower  altitudes,  (2)  the  fluorescence  of 
sunlight  at  2.7  urn,  and  to  a  lesser  extent,  2.0  and  1.6  urn  sunlight  in  the 
altitude  range  below  105  km,  and  (3)  the  transfer  of  excitation  from  oxygen 
singlet  0  to  vibration  via  a  quenching  collision,  0('D)  +  N2  is 
important  in  the  55  to  65  km  altitude  range.  Analysis  of  the  night  time 

4.3  ym  data  indicated  that  (1)  the  transfer  of  excitation  from  0H(V  >  0)  to 
N2(v«1)  is  required  to  model  the  data  in  the  75  to  90  km  altitude  region, 
and  (2)  a  second  mechanism  is  required  to  model  the  data  at  higher  altitudes. 
Direct  absorption  of  sunlight  and  a  "Chapman-like"  mechanism  were  found  to  be 
inadequate  for  this  purpose.  Our  initial  analysis  of  the  sunlit  2.7  ym  SPIRE 
data  stimulated  a  coordinated  effort  within  the  AFGL/DNA  community  to 
understand  this  aspect  of  the  SPIRE  data.  This  culminated  in  a  combined  and 
complimentary  effort  by  R.  Sharma,  T.  Degges  and  ourselves  in  arriving  at  an 
understanding  of  the  sunlit  2.7  ym  SPIRE  data.  The  major  mechanisms  were 
determined  by  this  triad  to  be  (1)  the  fluorescence  of  2.0  and  1.6  ym  sunlight 
by  CO2,  (2)  the  resonance  scattering  of  2.7  ym  sunlight  by  H20,  and  (3) 
absorption  and  reemission  of  2.7  ym  sunlight  by  C02. 

Quantitative  details  are  given  in  section  2.0  below,  previewed  by  a  detailed 
qualitative  summary  in  section  2.1. 

1.3  SUMMARY  OF  LABORATORY  EFFORT 

We  have  conducted  laboratory  studies  of  the  fluorescence  of  CO2  at  4.3  ym 
resulting  from  initial  excitation  of  the  101  band  due  to  absorption  of  laser 
light  near  2.7  ym.  By  utilizing  a  pulsed  laser  source  at  2.7  ym  we  have  been 
able  to  observe  the  resulting  fluorescence  at  4.3  ym  and  thereby  provide 
direct  laboratory  confirmation  of  this  mechanism  which  is  a  major  contributor 
to  the  earthlimb  radiance  from  C02  near  80  km  altitude. 


Our  report  discusses  several  experiments  designed  to  verify  the  branching 
ratio  for  radiation  at  4.3  pm  versus  reemission  at  2.7  pm.  Several 
experimental  difficulties  pertinent  to  this  measurement  are  discussed.  We 
also  describe  measurements  of  the  rate  of  decay  of  the  4.3  pm  emissions  for 

1  O  Ip 

CO2  and  CO2  in  the  presence  of  N2  and  the  interpretation  of 
these  results  as  they  relate  to  the  rate  of  vibrational  energy  transfer 
between  CO2  and  N2.  Additional  possibilities  for  enhanced  understanding 
of  the  mechanisms  related  to  CO2  emission  which  would  involve  measurements 
of  several  LWIR  bands  of  C0«  are  discussed. 


Section  2 

SPIRE  DATA  EVALUATION  AND  IMPLICATIONS 


2.1  DETAILED  QUALITATIVE  SUMMARY 

Analysis  of  the  15  ym  SPIRE  data  given  in  section  2.2  below  shows  that  our 
pre-SPIRE  non-LTE  model  for  the  C02  15  pm  earthlimb  was  deficient  in  two 
aspects.  First,  the  altitude  dependence  of  the  line  and  band  shape  needed  to 
be  accounted  for  in  our  radiative  transfer  model.  Second,  we  also  needed  to 
utilize  an  updated  rate  constant  for  the  quenching  of  C02(010)  by  atomic 
oxygen.  This  updated  rate  constant  is  required  to  be  much  faster  at 
mesopheric  temperatures  than  the  rate  recommended  by  Taylor  (1974)  on  the 
basis  of  the  extrapolation  from  high  temperature  laboratory  data  that  the  log 
of  the  rate  constant  is  proportional  to  T 

That  the  first  of  these  effects  might  be  important  was  first  suggested  at  a 
DNA  meeting  some  years  back  by  T.  Degges  in  discussion  of  effects  seen  in 
15  ym  zenith  radiance  data  that  had  been  obtained  in  a  rocket  borne  experiment 
on  3/22/73  (Rogers  et  al.,  1976).  In  subsequent  analysis  of  these  15  ym 
zenith  radiance  data  and  another  set  obtained  2/14/74  (Rogers  et  al.,  1977)  we 
incorporated  the  effects  of  altitude  dependence  in  the  line  and  band  shape. 

But  we  found  there  was  still  indication  that  some  further  improvements  were 
required  to  achieve  satisfactory  agreement  with  the  15  ym  zenith  radiance  data 
sets.  To  this  end  we  proposed  that  a  much  faster  rate  constant  at  mesopheric 
temperatures  for  the  quenching  of  C02(010)  by  atomic  oxygen  than  was 
recommended  by  Taylor  (1974)  on  the  basis  of  extrapolation  from  high 
temperature  laboratory  data  could  indeed  provide  the  key  to  achieving 
satisfactory  agreement  with  the  two  15  ym  zenith  radiance  data  sets. 

Recent  laboratory  measurements  (Simpson,  1980;  Allen  et  al.,  1980)  of 
quenching  of  002(010)  by  N2  have  shown  that  high  temperature  results  of 
this  rate  constant  could  not  be  extrapolated  to  low  mesopheric  temperatures  by 
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the  method  recommended  by  Taylor,  and  that  in  fact,  the  extrapolation  severely 
underpredicted  the  low  temperature  rate  constant,  as  we  had  suggested  might  be 
the  case  for  the  quenching  by  atomic  oxygen  as  based  on  our  analysis  of  the 
15  ym  zenith  radiance  data.  Based  on  the  Simpson  data,  Sharma  and  Nadile, 
(1980)  proposed  a  square  root  dependence  at  low  temperature  for  the  quenching 
by  atomic  oxygen,  normalized  to  achieve  a  reasonable  fit  to  the  SPIRE  15  ym 
data.  Sharma  developed  a  non-LTE  model  for  the  C02  15  ym  earthlimb  that 
utilized  both  the  altitude  dependent  line  and  band  shape  for  radiative 
transfer  for  the  C02  15  ym  radiation,  and  also  the  faster  rate  constants  at 
mesospheric  temperatures  for  quenching  by  N2,  02  and  most  importantly,  by 
atomic  oxygen.  He  utilized  this  model  to  achieve  a  much  better  understanding 
and  fit  to  the  SPIRE  15  ym  data  than  had  been  achieved  by  earlier  pre-SPIRE 
calculations.  In  the  work  we  report  here  in  section  2.2  we  provide 
independent  confirmation  for  the  results  obtained  by  Sharma. 

The  comparison  of  our  predictions  and  the  daytime  SPIRE  4.3  ym  data  that  we 
report  in  section  2.3  below  confirms  the  predictions  of  James  and  Kumer  (1973) 
that  the  fluorescence  of  2.7  ym  sunlight  by  C02  is  a  major  mechanism  for 
producing  the  daytime  4.3  ym  earthlimb  on  the  tangent  altitude  interval  from 
approximately  60  to  110  km.  The  role  of  fluorescence  of  2.7  ym  sunlight  by 
minor  isotopic  C02  species,  and  the  role  of  fluorescence  of  2.0  and  1.6  ym 
sunlight  by  C02  are  also  considered.  These  are  shown  to  be  significant,  but 
not  dominant  mechanisms  for  producing  the  daytime  4.3  ym  earthlimb.  The 
contribution  of  "weak"  C02  4.3  ym  bands  is  seen  to  be  approximately  as 
important  as  fluorescence  of  shorter  wavelength  sunlight  for  producing  the 
daytime  4.3  ym  earthlimb  for  tangent  altitudes  less  than  approximately 
110  km.  This  result  could  be  inferred  from  the  papers  of  Kumer  (1977a)  and 
Caledonia  et  al.  (1982).  For  tangent  altitudes  greater  than  approximately 
65  km  the  primary  weak  band  excitation  mechanism  is  absorption  of  4.3  ym 
sunlight.  This  is  also  the  primary  excitation  mechanism  for  the  strong  C02 
626  OOl-OOO  4.3  ym  band  in  that  altitude  region.  In  the  55  to  65  km  altitude 
region  the  excitation  of  N2(v*l)  by  transfer  of  energy  from  oxygen  singlet  D 
is  found  to  be  a  major  contributor  for  excitation  of  C02  4.3  ym  emission  as 
predicted  by  Kumer  (1977a).  It  was  found  that  a  better  fit  to  the  data  could 


be  achieved  in  the  altitude  region  less  than  80  km  by  utilizing  an  earthlimb 
radiance  calculation  that  accounts  for  altitude  dependence  in  line  and  band 
shape.  An  analytic  approximation  for  the  radiative  transfer  width  and 
transmission  functions  (defined  in  the  paper  by  Kumer  and  James,  1974)  derived 
from  work  done  by  Neuendorffer  (1982)  was  utilized  for  this  purpose.  The 
approximation  also  utilizes  the  mass  path  averaging  technique  to  approximate 
the  line  and  band  shape  parameters  for  a  given  earthlimb  viewing  line  of  sight 
as  suggested  by  Weinreb  (1982). 

In  general  we  slightly  overpredicted  the  daytime  4.3  urn  data  for  all  four 
sunlit  scans.  We  investigated  the  possibility  that  this  resulted  from  some 
uncertainty  in  the  knowledge  of  the  azimuth  of  the  SPIRE  scan  pattern  such 
that  the  pattern  was  rotated  further  than  thought  from  the  sun  thereby 
requiring  lower  solar  elevation  angles  for  our  calculation  than  were  actually 
prescribed  by  AFGL  personnel.  Careful  study  indicated  this  was  not  the  case. 
We  now  believe  our  general  overprediction  is  most  likely  the  result  of  a  minor 
inaccuracy  in  absolute  radiometric  calibration. 

The  analysis  of  the  night  time  4.3  ym  earthlimb  in  Section  2.4  supports  the 
theory  of  Kumer  et  al.  (1978)  that  the  principal  night  time  C02  4.3  ym 
excitation  mechanism  in  the  approximate  region  75  to  90  km  is  the  transfer  of 
vibration  from  the  OH  airglow  upper  states  0H(V  >  0)  to  N2(V=1)  followed  by 
subsequent  transfer  to  the  CO^vj)  states.  At  higher  altitude  an 

additional  source  of  excitation  is  required  to  model  the  data.  In  our  analysis 

we  showed  that  direct  absorption  of  sunlight  could  not  be  responsible  for  this 
additional  source.  We  also  showed  that  a  "Chapman-like"  airglow  related 
source  could  not  be  responsible  for  the  additional  excitation. 

In  section  2.5  we  discuss  our  work  in  modeling  the  C02  contribution  to  the 

sunlit  2.7  ym  earthlimb.  In  our  first  efforts  at  modeling  the  daytime  C02 

contribution  to  the  2.7  ym  earthlimb  we  showed  (Kumer,  1981)  that  absorption 
and  reemission  of  2.7  ym  sunlight  by  the  626  C02  specie  would  severely 
underpredict  the  SPIRE  data.  This  stimulated  the  community  to  reconsider  the 
problem.  Sharma  (1981)  showed  that  the  C02  contribution  would  be 


12 


9 


1 


considerably  enhanced  by  considering  fluorescence  of  2.0  and  1.6  ym  sunlight. 
Degges  (1981)  showed  there  was  a  very  significant  contribution  to  the  daytime 
2.7  ym  earthlimb  due  to  absorption  of  2.7  ym  sunlight  by  l-^O  followed  by 
subsequent  reemission.  We  showed  that  an  additional  enhancement  in  the  model 
could  be  obtained  by  considering  the  absorption  and  reemission  of  sunlight  at 
2.7  urn  by  the  minor  isotopic  CO2  species.  Sharma  (1981)  pulled  all  these 
results  together  and  showed  that  the  combination  of  these  effects  can  in  fact 
reasonably  account  for  the  SPIRE  daytime  2.7  ym  earthlimb  data. 

In  section  2.6  we  considered  a  simplified  implication  of  the  SPIRE  data. 
Looking  ahead  to  problems  of  interest  such  as  that  of  the  horizontal  structure 
in  limb  radiances  at  various  wavelengths  we  utilized  what  we  had  learned  in 
the  15  ym  region  to  generate  an  example  of  the  horizontal  structure  that  might 
be  caused  in  the  15  ym  region  by  a  gravity  wave.  For  the  simplified  gravity 
wave  model  we  utilized  a  cross  section  of  data  (Noxon,  1978;  Faire  et  al., 
1972;  Petersen  1979)  from  the  literature  on  the  perturbation  effects  of 
gravity  waves  on  the  ambient  temperature  profiles  and  the  simple  theory  of 
Hines  (1960)  to  relate  temperature  and  density  perturbations.  These  data 
suggested  that  a  model  with  a  temperature  perturbation  of  40  deg  at  95  km  and 
wavelengths  of  140  km  horizontally  and  20  km  vertically  is  not  unreasonable. 
This  example  gravity  wave  is  consistent  with  description  of  the  data,  and  the 
theory  of  the  phenomenon  presented  in  Chapter  5  in  the  report 
(AFGL-TR-81-0207,  Atmospheric  Infrared  Radiance  Variability,  edited  by 
Humphrey  et  al.,  1981),  which  is  devoted  to  structure  and  wave-like  phenomena 
in  the  upper  atmosphere.  For  this  model  we  perturbed  the  SPIRE  model 
atmosphere  for  two  cases  differing  by  *  in  phase  and  performed  15  ym  limb 
radiance  calculations  for  these  two  cases.  Significant  horizontal  radiance 
structure  was  noted  as  the  result  of  this  exercise. 

2.2  ANALYSIS  OF  SPIRE  15  ym  DATA 

The  discussion  in  section  2.2  focusses  on  the  analysis  of  the  final  version  of 
firm  computer  processed  spectrally  integrated  SPIRE  data  in  the  14  to  16  ym 
band  as  prepared  by  Sharma  and  Nadile,  1980.  The  final  version  of  the 
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Nadile  SPIRE  14  to  16  ym  data  and  a  representation  of  these  data  that  we  will 
use  in  subsequent  comparisons  with  our  calculations  is  shown  in  Fig.  2-1.  An 
overall  numerical  correction  factor  of  2,  discoverea  by  Dave  Green  in  1981, 
has  been  applied  to  the  data  shown  on  Fig.  2-1. 

2.2.1  Modifications  in  our  CO2  15  urn  Model  Required  to  Achieve  the  Best  Fit 
to  the  SPIRE  DATA 

A  major  point  of  interest  in  the  SPIRE  15  ym  data  is  the  prominent 
plateau-like  structure  which  appears  in  the  90-  to  110-km  region,  this  feature 
was  not  predicted  by  pre-SPIRE  non-LTE  models  of  the  15-ym  CO2  earthlimb. 

In  this  report  we  provide  evidence  for  upgrading  the  pre-SPIRE  model  by  two 
effects  (1)  inclusion  of  altitude  dependence  in  the  C02  rotational 
population  and  in  the  absorption  line  shape,  and  (2)  the  use  of  a  faster  rate 

O 

for  quenching  by  0(  P)  of  the  CO2  upper  state  01 '0  for  15-ym  radiation 
than  earlier  literature  recommended  estimates  (Taylor,  1974). 

The  principal  excitation  mechanisms  for  CO2  15-ym  emission  in  the  earth 
upper  atmosphere  include  absorption  of  earth  shine  and  excitation  by  thermal 
collisions  with  the  more  abundant  species  in  the  earth's  atmosphere,  i.e., 

N,  )  N- 

tf  /+  Ol'O  *  n  +  00°0  (2.1) 

U2  )  *  u2 

and 

0(3p)  +  Ol'O  *  0(3P)  +  00°0  (2.2) 

These  reactions  proceeding  to  the  right  result  in  the  quenching  component  due 
to  constituent  M  at  a  rate  Qijy)  =  Ki  [M]  where  M  may  be  02,  N2,  or 
O(^P),  and  where  i  correspondingly  is  1  or  2.  Proceeding  to  the  left 
excitation  by  thermal  collisions  is  given  by  Xim  =  2  k,Bk  [M]  where  Bk 
is  the  Boltzmann  factor  associated  with  the  local  atmospheric  kinetic 
temperature,  and  factor  2  accounts  for  the  increased  degeneracy  in  the  upper 
level  Ol'O  with  respect  to  the  lower  00°0.  For  our  discussion  it  is 
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critically  important  to  point  out  that  at  high  altitudes  (Z  ~  75)  where 
A/Q  >>  1,  i.e.,  where  radiation  rather  than  quenching  dominates  de-excitation 
of  the  Ol'O  state,  an  increase  in  the  quenching  rate  constants  actually 
results  in  an  increase  in  C02  15  pm  emission. 

Another  complicating  facet  of  the  problem  is  that  the  C02  15  urn  emission  is 
strongly  self-absorbed  along  the  limb  viewing  line  of  sight  up  to  120  km 
tangent  altitude.  Radiation  transport  that  accounts  for  multiple  absorption 
and  re-emission  of  CO 2  15  ym  photons  is  required.  We  will  review  our 
detailed  radiation  transport  formulation,  which  includes  the  temperature 
induced  altitude  variation  in  rotational  population  and  absorption  line 
shapes,  in  section  2.2.2  below. 

The  calculations  we  report  here  were  performed  with  a  model  atmosphere  for  the 
SPIRE  event  that  was  utilized  by  Sharma  and  Nadile  (1980). 

This  model  atmosphere  is  shown  on  Figs.  2-2  and  2-3.  The  Refs.  2,  3,  and  4 
that  are  cited  in  Table  2-1  are  given  in  Table  2-2.  Calculations  of  the  SPIRE 
14-  to  16-ym  earthlimb  above  65  km  were  performed  for  various  conditions  as 
listed  in  Table  2-1.  The  parameter  "key"  identifies  these  conditions  in 
subsequent  discussions.  The  parameter  Zb  is  the  lower  boundary  of  the 
calculations;  blackbody  radiation  upward  incident  at  Zb  is  characterized  by 
the  atmospheric  temperature  at  Zb.  The  rate  constant  k^  for  quenching  of 
Ol'O  by  or  02  and  k2  for  quenching  by  0  are  given  as  a  function  of 
temperature  T  by  the  following: 


Process  Reference 

Analytic  Approximations 

M  +  010  3  K,t 

M  +  010  4  kjj 

0+010  3  k?T 

0  +  010  2  k?s 

.  6.69  X  io:i°  *  l-si-om1'3) 

-  k1T  + 10"16  ^r 

«  2.32  x  10-9  x  exp  (-76.75/T1/3) 

-  k2T  +  io-14vr 
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LOG  10  OF  NUMBER  DENSITIES  (IN  cm 


Figure  2-3  Number  density  models  used  for  the  modeling  analysis  of 
the  SPIRE  data  (from  Sharma  and  Nadile, 6-27-1980) . 
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Table  2-1  CONDITIONS  USED  FOR  CALCULATIONS 


KEY  Zb 

(km) 

Rate  Constants 

M+010  0+010 

Absorption  Line  and 

Band  Shape(a) 

1 

65 

Ref.  4 

Ref.  2 

Variable  with  altitude 

2 

65 

Ref.  4 

Ref.  3 

Variable  with  altitude 

3 

65 

Ref.  3 

Ref.  3 

Variable  with  altitude 

3-LOLORP 

65 

Ref.  3 

Ref.  3 

Line  shapes  fixed  by  65-km 
temperature  and  65-km 
pressure  broadening 

3-HIL0RP 

65 

Ref.  3 

Ref.  3 

Line  shapes  fixed  by  65-km 
temperature  and  140-km 
pressure  broadening 

4-L0L0RP 

65 

Ref.  4 

Ref.  2 

As  in  3-LOLORP 

4-HILORP 

65 

Ref.  4 

Ref.  2 

As  in  3-HIL0RP 

5 

40 

Ref.  4 

Ref.  2 

Variable  with  altitude 

6-HILORP 

40 

Ref.  4 

Ref.  2 

Line  shapes  fixed  by  40-km 
temperature  and  140-km 
pressure  broadening 

(a)  Distribution  of  rotational  populations. 


Table  2-2  REFERENCES  FOR  RATE  CONSTANTS  UTILIZED  IN  THE  CO£  15  vm  MODEL 


2.  R.  Sharma  and  R.  Nadile,  "Carbon  Dioxide  (v2)  Radiance  Results  Using 

a  New  Nonequi librium  Model,"  paper  presented  at  DNA-AFGL  Nuclear  Weapon 
Effects  Chemistry  Conf.,  AFGL,  Bedford,  MA,  26  and  27  June  1980. 

3.  R.  L.  Taylor,  "Energy  Transfer  Processes  in  the  Stratosphere,"  Can.  J. 
Chem. ,  Vol.  52,  1974,  p.  1436. 

4.  C.  J.  S.  M.  Simpson,  private  communications,  1980. 
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The  analytic  approximations  to  and  k2$  were  presented  by  Sharma  at  the 
1980  DNA  Chemistry  Meeting.  The  analytic  forms  for  and  K2$  shown  here 
are  preliminary  expressions  derived  by  Sharma.  The  form  for  was 
selected  by  Sharma  to  have  the  same  low  temperature  functional  form  as  K1S 
and  to  achieve  a  reasonable  fit  to  the  SPIRE  data.  We  stress  that  the  form 
that  we  show  for  is  a  preliminary  version  and  that  Sharma  has  probably 

updated  this  version  in  his  more  recent  efforts  to  achieve  a  best  fit  to  the 
SPIRE  data. 

Calculations  in  which  variation  in  the  C02  band  and  line  absorption  shapes 
as  a  function  of  altitude  varying  temperature  and  pressure  are  designated  by 
key  =  1,  2,  3,  and  5.  For  comparison,  some  calculations  performed  with  band 
and  line  shapes  constant  throughout  the  atmosphere  are  characterized  by 
temperature  T  taken  to  be  T(Zb)  and  by  pressure  broadening  as  given  by 
Table  2-1. 

Case  key  =  3-HILORP  represents  the  pre-SPIRE  model,  which  uses  the  old  rate 
constants  and  non-varying  line  and  band  shapes;  results  of  the  calculation  are 
compared  to  the  data  representation  in  Fig.  2-4.  It  can  be  seen  that  the 
plateau-like  feature  present  in  the  data  near  95  to  110  km  is  not  at  all 
present  in  the  calculation.  The  case  key  =  2,  which  improves  on  the  old  model 
in  that  the  recent  Simpson  data  for  reaction  (1)  is  utilized  and,  in  that  the 
altitude  variation  in  line  and  band  shape  is  accounted  for,  is  shown  in 
Fig.  2-5.  Here  we  see  noted  improvement  between  the  predictions  and  data  as 
compared  to  the  pre-SPIRE  results  previously  shown  on  Fig.  2-4.  The  best 
agreement  however  is  achieved  for  the  case  key  =  1  and  is  shown  on  Fig.  2-6. 
This  case  further  includes  a  postulated  low  temperature  enhancement  in  0(Jp) 
rate  constant  k ^  for  quenching  Ol'O  as  discussed  above.  We  consider  the 
case  key  =  1  as  providing  the  best  fit  since  details  of  the  shape  of  the  limb 
radiance  profile  are  best  reproduced  by  this  case.  A  somewhat  higher  quality 
fit  of  the  data  could  be  achieved  in  this  case  by  translating  the  SPIRE  data 
set  down  by  -  2  to  4  km.  This  is  consistent  with  the  stated  experimental 
accuracy  of  absolute  altitude  determination  for  the  SPIRE  data. 
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KEY  =  3.  HILORP 


TANGENT  ALTITUDE  (km) 


Figure  2-4  Calculation  based  on  the  pre-SPIRE  model  which  uses  the  old 
rate  coefficients  and  altitude  independent  line  and  band  shapes  is 
compared  with  the  SPIRE  data,  key* 3 -HILORP. 


KEY  =  2, 


TANGENT  ALTITUDE  (km) 


Figure  2-5  Calculation  improved  on  pre-SPIRE  model  in  that  altitude 
dependence  of  line  and  band  shape  is  accounted  for, the  old  rate 
coefficients  are  utilized  however,  is  compared  with  the  SPIRE  data, 
key«2. 


KEY  =  1, 


TANCENT  ALTITUDE  (km) 


Figure  2-6  Final  version  of  the  calculation  is  compared  with  the  SPIRE 
data  ,  key*l.  This  version  is  improved  over  the  pre-SPIRE  model  in 
that  (1)  altitude  dependence  of  the  line  and  band  shapes  is  accounted 
for rand  (2)  the  newer  values  for  the  rate  constants  as  proposed  by 
Sharma  (discussed  in  the  text)  are  utilized. 


In  this,  our  final  report,  we  present  results  for  the  cases  key  -  1,  2,  and 
3-HILORP  since  these  consistute  a  useful  illustration  of  the  improvements  we 
have  implemented  in  our  15  ym  model.  Exhaustively  detailed  results  and 
further  discussion  of  all  the  cases  included  in  Table  2-1  have  been  presented 
in  a  previous  contractual  progress  report.  These  results  and  the  further 
discussion  are  available  on  request  for  the  interested  reader. 

2.2.2  Details  of  the  Improved  C02  15  um  Radiation  Transport  Model 

The  ability  to  account  for  the  temperature  and  pressure  altitude  dependent 
variation  in  the  CO^  010-000  line  and  band  shape  is  one  of  the  two  major 
improvements  in  the  models  that  are  required  to  achieve  an  improved  fit  to  the 
data.  In  this  sub-section  we  discuss  the  implementation  of  the  effect  into 
the  model  in  more  detail. 

The  radiation  transport  problem  breaks  into  two  distinct  steps:  first,  the 
non-LTE  calculation  of  the  altitude  dependent  population  of  upper  state  of  the 
band,  and  second,  the  calculation  of  the  observable  radiance  that  results  from 
this  non-LTE  solution.  In  step  1,  multiple  self  absorption  and  emission  of 
photons  must  be  accounted  for  in  general.  In  step  2,  emission  and  self 
absorption  of  photons  along  the  sensor  line  of  sight  must  be  accounted  for. 

Calculation  of  the  Non-LTE  Solution.  On  Fig.  2-7  we  show  the  basic  geometry 
and  transport  equation  to  determine  S  the  excitation  rate  of  C02  01  * 0  per 
unit  a  [CO^]  at  altitude  Z.  The  discretized  matrix  version  of  this  integral 
radiative  transport  equation  Cj-  Sj  =  SQi  which  actually  must  be 
solved  on  a  digital  computer  in  order  to  solve  for  Si  (i.e.,  the  value  of  S 
at  altitude  Z^)  is  also  shown  on  Fig.  2-7. 

The  physics  involved  in  the  calculation  of  the  radiative  transport  matrix 
elements  are  shown  in  Fig.  2-8.  The  nested  loops  in  the  computer  program  that 
are  required  for  the  calculation  of  the  matrix  elements  H i j  and  for  the 
array  Tj  is  also  indicated  on  Fig.  2-8. 
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BASIC  TRANSPORT  EQUATION 


Figure  2-7  Geometry  and  formulation  of  the  basic  non-LTE 
radiation  transport  equations 


OMPUTE  H.. 


Figure  2-8  Computation  of  the  radiation  transport  equation 

kernal  function  matrix  elements  including  temperature 
and  pressure  variation  in  line  and  band  absorption  profiles 


Note  that  line  shape  altitude  dependence  is  accounted  for  by  the  altitude 
dependent  parameters  ai  and  Gj(Z.j).  The  quantity  bi  is  the 
ratio  of  Doppler  width  vq(Z|3)  at  reference  altitude  Zb  to  the  width 
vD(Zi)  at  the  altitude  Z-j.  Here  vq  is  the  1/2  e-fold  Doppler  width  as 
used  commonly  in  the  astrophysical  literature.  The  altitude  dependent 
parameter  a^  gives  the  ratio  of  the  1/2  maximum  1/2  width  pressure  broadened 
Lorentz  line  parameter  P(Z.)  Y  (T[zi 3/296°)_1/2  to  the  Doppler  width 
vD(Zi).  Here  P(Zi)  is  the  atmospheric  pressure  at  altitude  Z^  in 
units  atm  at  altitude  Z^,  Y  ~  0.07  cnf^/atm  is  the  pressure  broadened 
Lorentz  parameter  for  C02  by  air  at  STP  and  the  factor  (T(Zi )/296°)"1/2 
accounts  for  temperature  variation  in  y.  The  Voigt  line  profile  tf(a.j,  b^ 
x)  then,  with  x  =  (v  -  v0)/vD(Zb),  gives  the  altitude  dependence  of  the 
line  shape  for  a  line  centered  at  vQ.  The  altitude  dependence  of  the  band 
shape  is  given  by  Gj(Z.j),  the  fractional  rotation  population  of  C02  in 
the  level  J.  The  quantities  GjtTtZ^]  and  v^T^)]  are  explicitly 
dependent  on  atmospheric  kinetic  (and  rotational)  temperature  T(Z.j),  and 
derive  their  altitude  dependence  as  the  result  of  altitude  dependence  of  the 
temperature. 

Also  note  the  exponential  function  (SbT  jJ(x))which  occurs  in  the 
expression  for  H.^;  this  function  relates  emission  and  absorption  between 
layer  pairs  i  and  j  in  a  plane  parallel  atmosphere.  This  function  has  been 
defined  and  discussed  in  application  to  radiative  transport  in  plane  parallel 
atmospheres  going  back  to  the  days  to  Milne.  The  property  that  E^  (r)  is 
proportional  to  -In  (r)  forT  —  0  requires  the  use  of  the  somewhat  complicated 
expression  for  as  given  on  Fig.  2-7.  This  expression  exploits  the  fact 
that/^dT-'  E^  (r‘ )  =  E2  (o)  -  E2  (t)  where  E2(o)  =  0.5  and 
E2(t)  decreases  monotonical ly  with  increasing  t.  This  property  is  exploited 
by  utilization  of  the  escape  function  Ei  as  defined  on  Fig.  2-9.  Note  the 
alternate  interpretation  that  Ei  =  L-  H..  i„  the  limit  that  the  j  mesh 
goes  toco.  Thus  on  a  finite  altitude  mesh  we  can  easily  compute  well  behaved 
matrix  elements  G^  *  -IH^  Wi  ffj  for  i  *  j  ana  Gi  i  =  1  -  (1  - 
+Ii*j  Hij  wjfij*  The  <luant’ty  Qj  is  the  fraction  of 
de-excitation  of  the  upper  level  at  Zi  due  to  radiation.  Total 
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PROBABILITY  OF  ESCAPE  FROM  LEVEL  i 


Figure  2-9  Computation  of  escape  function  and  transmission 
function  matrix  elements 


de-excitation  is  due  to  the  sum  of  the  radiative  and  quenching  processes. 

This  is  consistent  with  the  definition  of  as  proportional  to  the  total 
excitation  rate  at  Zi  in  our  formulation  of  the  problem.  The  quantities 
W.  are  weighting  factors  required  to  perform  the  radiative  transport 
integrals  over  altitude.  More  detail  on  the  choice  of  the  W-  is  given  in 

v 

the  paper  by  Kumer  and  James  (1974).  Note  here  that  v.^.  =  y*!°  0 
r  J  J  ^b 

[C02]  dz  =  a  where  Gj  Sg  o  is  ^  the  zenith  optical  depth  at 

altitude  Zb  at  line  center  for  a  line  in  branch  B  (where  B  can  be  p,  q  or  r) 

which  originates  from  ground  rotational  level  J.  Detailed  definitions  and 

properties  of  the  Ep  (r)  can  be  found  in  standard  texts  such  as  Abramowitz 

and  Stegun  (Handbook  of  Mathematical  Functions.  U.  S.  Dept,  of  Commerce,  U. 

S.  Printing  Office,  Washington,  D.  C.  1966). 

The  element  then  is  formed  by  an  integral  over  all  lines  in  the  band, 
and  over  all  frequencies  in  each  line,  it  gives  the  contribution  to  Sj  due 
to  absorption  at  level  Zj  of  photons  that  are  emitted  from  level  Z^.  The 
calculation  assumes  that  the  lines  in  the  band  are  non  over  lapping,  that  gj 
applies  for  both  the  lower  and  upper  states  of  the  band,  and  that  i>  is  the 
line  profile  for  both  absorption  and  emission  at  1-. 

J 

To  calculate  the  solution  to  the  discretized  radiative  transport 
equation,  equation  G-jj  Sj  =  Soi  then  requires  a  computation  of  the 
matrix  elements  as  described  here  and  on  Figs.  2-7  and  2-8,  and  a  computation 
of  the  inhomogeneous  term  SQi.  This  followed  by  a  routine  calculation  of 
the  inverse  matrix  G~1  allows  for  calculation  of  the  via  S,  = 

I J  II 

6ij  Soj- 

Given  the  it  is  easy  to  use  standard  matrix  inversion  routines  that  have 
been  developed  over  the  years  by  many  different  groups  for  matrix  inversion  to 
obtain  the  elements  G'jj. 

The  inhomogeneous  elements  SQi  are  the  result  of  absorption  of  earthshine 
which  crosses  the  lower  boundary  Zb,  and  a  second  mechanism,  collisional 


29 


excitation.  This  second  mechanism  was  discussed  in  considerable  detail  in 
section  2.2.1. 

The  details  of  calculating  the  contribution  to  SQi  due  to  the  absorption  of 
earthshine  are  given  in  our  earlier  papers  (Kumer  and  James,  1974;  Kumer  1977a 
and  1977b).  In  the  present  application  we  have  upgraded  this  calculation  to 
include  the  altitude  dependence  of  the  line  shape  i>  (a^  b^  x)  and  band 
shape  Gj  (Z.,). 

The  Radiance  Calculation.  Once  the  non-LTE  solution  S.  has  been  computed 
the  next  step  is  to  calculate  the  radiance,  as  might  be  observed  by  a  rocket 
borne  sensor  in  either  a  zenith  or  an  earthlimb  geometry,  as  is  shown  in 
Figs.  2-9  and  2-10  respectively.  The  formulation  shown  utilizes  altitude 
dependent  line  and  band  shapes. 

2.2.3  Further  Improvements  Required  in  Our  15  pm  C0->  Model 

To  date  we  have  implemented  the  altitude  dependent  line  and  band  shape 
formulation,  as  described  in  section  2.2.3  above,  for  just  the  C02  15  pm 
fundamental  band  0r0*00°0.  In  our  modeling  effort,  the  results  of  which 
are  shown  on  Figs.  2-4,  2-5,  and  2-6,  we  found  that  hot  bands,  such  as 
0220  Ol'O,  also  contribute  significantly  to  15  pm  C02  radiance.  For  these 
bands  however  we  remain  limited  to  our  older  formulation  which  does  not 
account  for  altitude  dependence  in  the  line  and  band  shapes. 

In  order  to  perform  detailed  spectral  modeling  of  the  SPIRE  15  pm  data,  and  of 
data  we  anticipate  will  be  obtained  in  the  upcoming  experiments  such  as  CIRRIS 
and  SPIRIT,  we  will  need  to  upgrade  our  code  to  account  for  this  effect  in  all 
the  bands.  The  code  also  needs  to  be  upgraded  so  that  all  isotopes  can  be 
dealt  with  in  a  single  run  rather  than  a  number  of  separate  runs. 

At  lower  altitudes  the  radiance  calculation  needs  to  be  modified  to  account 
for  line  overlap.  We  will  show  explicit  examples  of  this  effect  and  how  it 
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LIMB  OPTICAL  DEPTHS 


impacts  analysis  of  4.3  pm  SPIRE  data  in  section  2.3.  The  effect  will  be 
important  in  low  altitude  15  pm  data  as  well,  particularly  in  the  Q-branch 
region. 


2.3  ANALYSIS  OF  SUNLIT  SPIRE  4.3  pm  DATA 
2.3.1  Data  Processing  and  Comparison 

To  assist  us  in  this  analysis  Nadile  and  Grieder  provided  us  with  spectrally 
integrated  SPIRE  data  on  the  following  spectral  pass  bands. 


Pass  Band  (pm) 

Spectrometer 

■  m 

-  4.03 

HS-2 

989 

-  4.45 

HS-2 

4.45 

-  4.48 

HS-2 

3.90 

-  4.10 

NS-2 

4.10 

-  4.50 

NS-2 

A  comparison  of  our  prediction  for  spectrally  integrated  4.3  pm  earthlimb 
radiance  with  the  spectrally  integrated  SPIRE  NS-2  data  for  sunlit  scans  9 
through  12  are  shown  on  Figs.  2-11  through  2-14.  The  data  (curve  1  on 
Figs.  2-11  through  2-14)  are  given  by 

D  =  (0.4  pm)  (R4 el  _  415  -  R3>9  _  4>1) 

where  R4i]_4<5  is  the  average  spectral  radiance  from  4.1  to  4.5  pm  provided 
to  us  by  Nadile  and  Grieder  and  R3>g_4>i  is  the  average  spectral  radiance 
from  3.9  to  4.1  pm.  The  R3>g_4ti  is  subtracted  from  R4,i_4,5  in  order  to 
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.TANGENT  ALTITUDE  IN  KM 


Figure  2—11  The  date  (curve  1)  for  scan  9  are  compared  with  our 
predictions (curve  2)  .Curve  3  is  the  component  of  the  predictions  du€ 
to  emission  by  the  626  CO2  001-000  bandfcurve  4  is  the  component  due 
to  the  "weak"  bands  and  curve  5  the  solar  flourescence  component. 


provide  a  first  order  correction  for  DC  offset  in  the  spectral  data.  These 
predictions  for  scans  9  through  12  apply  for  solar  elevation  data  that  were 
also  provided  to  us  by  Nadile  and  Grieder.  The  solar  elevation  data  were  a 
subset  of  comprehensive  SPIRE  attitude  and  trajectory  data  that  were  provided 
by  Nadile  and  Grieder.  The  solar  elevation  angle  at  the  tangent  point  of  the 
sensor  line  of  sight  on  the  40  to  140  km  tangent  altitude  range  are  listed  for 
scans  9  through  12  below: 


Scan  No. 

Range  of  Solar  Elevation  (deg) 

9 

13.2  to  9.6 

10 

11.6  to  7.3 

11 

6.8  to  4.2 

12 

10.2  to  0.3 

2.3.2  Spectral  Composition  of  the  Sunlit  C0?  4.3  pm  Earthlimb 

Our  predictions  for  the  sunlit  C02  4.3  ym  earthlimb  are  labelled  by  curve  2 
on  Fig.  2-11  through  2-14,  and  are  broken  into  components  that  are  curves  3 
through  5.  Curve  2  is  the  sum  of  these  components.  Some  of  the  bands  that  we 
have  included  in  our  predictions  of  the  sunlit  4.3  ym  earthlimb  are  listed  in 
Table  2-3.  Curve  3  is  the  contribution  from  the  fundamental  band  of  the  major 
CO2  isotope,  i.e.,  the  626  001-000  band.  Curve  4  is  the  component  that 
includes  the  "weak"  bands  no.  2  through  8  in  Table  2-3.  The  principal 
excitation  mechanism  for  the  components  labelled  by  curves  3  and  4  is  the 
resonance  scattering  of  sunlight,  i.e.,  the  absorption  of  sunlight  at  4.3  ym 
by  CC>2  followed  by  subsequent  reemission  at  4.3  ym.  A  number  of  other 
mechanisms  contribute  to  the  day  time  4.3  ym  earthlimb  as  was  discussed  by 
Kumer  and  James  (1974).  These  are  listed  in  Table  2-4  and  will  be  discussed 
below  in  more  detail  in  the  perspective  of  their  importance  to  understanding 
the  SPIRE  data.  Curve  5  is  the  solar  fluorescence  contribution  which  is  due 
to  absorption  of  sunlight  by  CO2  at  wavelengths  shorter  than  4.3  ym  followed 
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Table  2-3  THE  C02  4.3  ym  BANDS 


j 

co2  (jv3)  (j) 

Transition 

Isotope  Comment 

1 

00°  1 

00°0 

626 

»  Strong  Band 

2 

00°  1 

00°0 

636  \ 

3 

00°  1 

00°0 

628 

4 

00°  1 

00°0 

627 

5 

oih 

01*0 

626 

>—  Weak  Bands 

6 

oih 

01*0 

636  | 

7 

oih 

01*0 

628 

8 

011! 

01*0 

627 

9 

f-H 

C\i 

00 

o 

0220 

626  > 

10 

02°  1 

02°0 

626  1 

j— Fluorescent  Bands 

11 

on 

io°o 

626  J 

by  reemission  at  4.3  ym.  The  contributions  to  curve  5  include  bands  10  and  11 
in  Table  2-3.  These  are  the  contribution  from  the  626  isotope  that  result 
from  absorption  of  sunlight  at  2.7  ym  followed  by  reemission  at  4.3  ym.  Curve 
5  also  includes  contributions  from  the  analogous  minor  isotopic  bands  that  are 
not  listed  in  Table  2-3.  Other  bands  which  have  been  included  in  the 
calculation  labelled  curve  5  are  those  C02  4.3  ym  fluorescence  bands  that 
result  from  absorption  of  sunlight  at  wavelengths  2.0  and  1.6  ym.  The 
relative  importance  of  these  bands  will  be  discussed  in  more  detail  below. 

Verification  that  solar  fluorescence  of  2.7  ym  sunlight  is  a  major  source  of 
4.3  ym  earthlimb  emission  was  a  driving  factor  for  performing  the  SPIRE 
experiment.  The  importance  of  this  mechanism  was  predicted  by  James  and  Kumer 
(1973).  It  is  seen  that  the  fluorescence  is  in  fact  the  major  contribution 
near  80  km  as  was  predicted. 
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Table  2-4  MECHANISMS  FOR  DAYTIME  EXCITATION  OF  C02  v3  IN  THE 
UPPER  ATMOSPHERE 


Item 

Mechanism 

1 

Thermal  collisions  with  atmospheric  species  other  than  N2.* 

2 

N+N0->N2(v>0)  +  0 

3 

0(lD)+N2->0(3p)+N2(V>0) 

4 

Thermal  collisions  with  N2 

5 

Solar  photo  electron  collisional  excitation 

6 

C02( j)+hv(4.3  ym)->C02( JV3) ,  i.e.,  absorption  of  4.3  urn 
sunlight  in  the  C02  "weak"  bands 

7 

626(000)+hv(4.3  pm)->626(001) ,  i.e.,  absorption  of  4.3  p m  sunlight 
in  the  fundamental  band  of  the  major  isotopic  specie  of  C02 

8 

Absorption  of  4.3  pm  earthshine  by  atmospheric  C02 

9 

Absorption  of  sunlight  by  C02  at  wavelengths  shorter  than  4.3 
pm,  mainly  at  2.7  pm  but  absorption  at  wavelengths  2.0  and  1.6  pm 
have  also  been  included,  to  produce  the  combination  states 

C02(Jv3),  next  the  J  portion  of  the  combination  state  is 
passed  to  C02  via  a  resonance  collisional  transfer,  i.e. 

C02(Jv3)+C02->C02(J)+C02(001) 

,  or  C02(001)  is  obtained  by  the  circuitous  process; 

C02( Jv3) +N2->C02( J)+N2( V=1 ) 

N2(v=1)+C02->N2+C02(001) 

,  in  each  case  the  end  product  of  interest  is  C02(001). 

♦The  mechanisms  above,  exclusive  of  mechanism  6,  have  been  discussed  in 
detail  by  Kumer  and  James  (1974);  mechanism  6  has  been  discussed 
in  detail  in  papers  by  Kumer  (1977ab). 
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It  is  also  seen  that  the  "weak"  bands  contribute  approximately  the  same  amount 
to  the  limb  radiance  near  80  km  as  does  the  fluorescence  component.  Both  the 
"weak"  band  and  the  fluorescence  emission  dominate  the  contribution  due  to  the 
626  001-000  band  for  altitudes  less  than  100  km.  The  "weak"  band  contribution 
to  the  sunlit  CO^  4.3  urn  earthlimb  emission  was  not  addressed  in  the 
original  James  and  Kumer  (1973)  paper,  but  its  importance  could  be  inferred 
from  the  later  papers  by  Kumer  (1977a  and  1977b).  Analogously,  Caledonia, 
Green,  and  Murphy  (1982)  point  out  that  the  weak  bands  are  the  dominant 
contributors  to  CO^  4.3  pm  emission  that  is  observed  in  laboratory 
experiments  that  involve  heavily  self  absorbed  CO2  mass  path  lengths.  Here, 
our  analysis  of  the  SPIRE  data  indicates  that  the  4.3  pm  CO^  mass  path  for 
earthlimb  geometry  is  large  enough  to  produce  significant  effects  for 
altitudes  less  than  or  equal  to  approximately  110  km.  The  pertinent  effect  is 
that  the  "weak"  bands  make  up  almost  1/2  of  the  4.3  pm  daytime  earthlimb 
emission  on  the  tangent  altitude  interval  from  approximately  55  to  95  km. 

2.3.3  General  Discussion  of  the  Comparison  of  the  Data  and  the  Predictions 

Here  we  note  that  although  the  data  vary  by  more  than  2  orders  of  magnitude  on 
the  range  40  to  140  km  tangent  altitude,  our  predictions  are  never 
significantly  more  than  a  factor  of  2  in  disagreement,  and  for  the  most  part 
agree  better  than  that  with  the  data.  There  are  several  cases  in  the  overall 
comparison  of  the  predictions  and  data  for  the  four  scans  where  disagreement 
is  more  obviously  marked,  and  we  will  discuss  these  cases  here  in  general,  and 
then  in  subsequent  discussion  in  more  detail. 

First  we  note  relatively  large  disagreement  between  prediction  and  data  in  the 
high  altitude  end  of  scans  10  and  11.  Inspection  of  the  predictions  for  all 
four  scans  and  of  the  data  for  the  presumably  better  behaved  scans  (9  and  12) 
would  argue  that  the  discrepancy  can  be  attributed  to  noise,  or  perhaps  an 
effect  of  the  variation  of  the  DC  level,  and  therefore,  should  be  discounted. 
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Secondly,  we  note  a  sharp  disagreement  between  prediction  and  data  at  lower 
tangent  altitudes  beginning  at  54  km  and  becoming  very  noticeable  at  40  km  in 
all  cases.  We  are  reasonably  sure  that  this  results  from  a  breakdown  in  our 
assumption  of  non-overlapping  lines  in  the  calculation  of  earthlimb  radiance 
for  tangent  altitudes  less  than  -  55  km.  We  discuss  this  point,  and  a  method 
to  remedy  the  situation  in  detail  in  section  2.3.7  below. 

Third,  there  appears  to  be  a  general  tendency  for  over  prediction  of  varying 
significance  on  the  entire  tangent  altitude  range  from  55  to  140  km  for  the 
four  scans.  This  apparent  discrepancy  is  not  particularly  serious  (typically 
much  less  than  a  factor  2)  and  here  we  briefly  consider  several  alternative 
possibilities  to  improve  on  this  situation;  (a)  the  generally  mild 
overprediction  on  the  approximate  range  ~  80  to  140  km  would  be  reduced  if  we 
would  use  the  Hays  and  Olivero  (1970)  C02  mixing  ratio  model,  which  breaks 
away  from  and  becomes  less  than  the  constant  altitude  value  of  3.2  x  1CT4  at 
80  km,  rather  than  at  100  km,  which  is  the  model  we  have  employed  for 
generating  our  predictions  on  Figs.  2-11  through  2-14;  (b)  next,  there  is  a 
possibility  that  the  solar  flux  is  less  than  the  value  that  we  used  in  the 
calculation,  i.e.,  there  may  be  solar  absorption  lines  coincident  with  the 
C02  atmospheric  absorption  lines  near  4.3  ym;  we  have  considered  this 
possibility  to  the  extent  that  we  have  looked  at  papers  such  as  thatof  Mount 
and  Linsky  (1975)  which  indicate  that  it  is  very  doubtful  that  there  is  enough 
C02  in  the  solar  atmosphere  to  cause  the  required  absorption  in  the  lines  of 
the  fluorescent  and  weak  bands  that  dominate  the  4.3  ym  earthlimb  from  ~  50  up 
to  -  105  km,  and  have  concluded  that  this  explanation  is  very  unlikely;  (c) 
next,  we  have  also  considered  the  possibility  that  the  entire  scan  pattern  of 
the  SPIRE  experiment  was  rotated  further  away  from  the  sun  than  was  determined 
by  AF6L  personnel  on  the  basis  of  star  sensor  and  other  data.  Here  we 
consider  item  (c)  in  summary. 

We  discuss  the  issue  of  possible  azimuthal  misregistration  in  detail  in 
section  2.3.8  below.  Attitude  sensor  systems  on  the  SPIRE  did  experience  some 
problems.  Vertical  positioning  was  solved  to  better  than  5  km  accuracy  on  all 
the  scans  by  the  utilization  of  prominent  and  well  known  features  in  spectral 


earthlimb  emission.  The  comparison  between  data  and  predictions  on  Figs.  2-11 
through  2-14  indicates  generally  that  vertical  registration  of  the  SPIRE  data 
was  indeed  better  than  5  km  in  accuracy.  The  general  overprediction  seen  in 
Figs.  2-11  through  2-14  suggests  the  average  in  whatever  vertical 
misregistration  that  there  may  have  been  should  shift  the  data  upwards 
relative  to  the  predictions. 

The  problem  is  somewhat  more  difficult  for  the  azimuthal  case.  We  looked  into 
this  possibility  in  some  detail,  and  report  on  it  in  section  2.3.8  below. 

Here  we  summarize  that  our  study  was  not  absolutely  conclusive  but  seemed  to 
indicate  that  a  major  error  in  azimuthal  registration  of  the  SPIRE  scan 
pattern  is  not  likely.  Again,  the  overprediction  noted  in  Figs.  2-11  through 
2-14  suggests  that  if  there  is  any  average  misregistration  in  the  SPIRE 
azimuth,  it  should  be  away  from  the  sun  to  reduce  the  predictions  relative  to 
the  data. 

Finally  we  believe  that  the  most  likely  explanation  for  the  general 
overprediction  might  be  the  simple  result  of  some  inaccuracy  in  the 
radiometric  calibration  of  the  SPIRE  instrument,  and  that  the  small  scan  to 
scan  variations  in  the  qualitative  comparison  between  the  data  and  the 
predictions  might  be  the  effect  of  some  minor  scan  to  scan  misregistration  in 
both  the  vertical  and  in  the  azimuthal  registration  of  the  altitude  of  the 
SPIRE. 


2.3.4  Mechanisms  for  C02  4.3  pm  Fluorescence 

In  discussion  above  it  has  been  pointed  out  that  4.3  pm  fluorescence  of 
sunlight  is  a  major  mechanism  for  the  4.3  pm  earthlimb  on  the  ~  60  to  ~  100  km 
tangent  altitude  range.  Here  we  discuss  that  mechanism  in  a  little  more 
detail.  The  James  and  Kumer  (1973)  work  was  limited  to  just  the  C02  626 
fluorescence  originating  from  absorption  of  2.7  pm  sunlight.  Recently  Sharma 
and  Nadile  (1981)  have  shown  that  fluorescence  of  ~  2.0  pm  and  -  1.6  pm 
sunlight  is  important  for  explaining  the  2.7  pm  SPIRE  data.  This  suggested 
that  we  could  improve  on  our  previous  work  by  quantitatively  accounting  for 
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this  effect  at  4.3  pm.  Also  it  has  been  seen  above,  and  it  also  has  recently 
been  recognized  by  a  number  of  workers  in  the  community,  that  contributions 
from  minor  isotopic  species  such  as  C02  636,  628,  and  627  can  be  very 
important  in  the  formation  of  C02  earthlimb  emission. 

In  the  work  we  report  here  we  have  accounted  for  these  effects  on  the  4.3  pm 
earthlimb  fluorescence.  On  Fig.  2-15  we  show  for  scan  9  the  contributions  to 
4.3  pm  earthlimb  fluorescence  from  these  sources.  Curve  1  on  Fig.  2-15  is  the 
contribution  from  fluorescence  of  2.7  pm  sunlight  by  C02  626,  i.e.,  the 
mechanism  considered  previously  by  James  and  Kumer.  Curves  2  through  4 
illustrate  the  contribution  of  the  C02  isotopic  species  636,  628,  and  627 
respectively.  Curves  5  and  6  show  the  contributions  due  to  C02  fluorescence 
of  2.0  and  1.6  pm  sunlight  respectively. 

It  can  be  seen  that  these  additional  sources  of  C02  4.3  pm  fluorescence  can 
not  be  completely  neglected.  Curve  5  on  Figs.  2-11  through  2-14  above  is  the 
sum  of  the  six  components  shown  on  Fig.  2-15.  Although  these  additional 
sources  of  fluorescence  as  defined  in  this  section  are  not  particularly 
important  for  the  4,3  pm  earthlimb,  the  analysis  by  Sharma  shows  that  they 
should  indeed  be  of  major  significance  for  the  sunlit  2.7  pm  earthlimb. 

2.3.5  Mechanisms  for  the  Sunlit  Earthlimb  Emission  at  4.3  pm 

The  mechanisms  for  excitation  of  C02  (Ov-j)  are  shown  in  Table  2-4.  The 
contribution  for  scan  9  of  each  mechanism  to  the  quantity  S ( 6  N)  (which  is 
proportional  to  the  excitation  rate  of  the  upper  state  626  001  of  the 
fundamental  C02  4.3  pm  band  of  the  major  isotope)  is  shown  on  Fig.  2-16.  We 
see  that  absorption  of  sunlight  (in  one  way  or  another)  is  the  dominant 
mechanism  all  the  way  down  to  approximately  55  to  ~  65  km  where  the  oxygen 
singlet  D  mechanism  begins  to  become  dominant.  The  contribution  by  the  0 
(*D)  mechanism  to  the  626  001-000  band  limb  radiance  is  shown  on  Fig.  2-17. 

It  is  seen  that  evidence  for  any  serious  discrepancy  between  our  model  for 
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CONTRIBUTION  TO  001-000  EMISSION  BY  MECHANISM 


th?m,ah2Q1LT5!K?°noribi!fci0I!S  from  each  of  the  mechanisms  numbered  1 
through  9  on  table  2-4  to  the  quantity  ,  626  001-000  4.3um 

emission/ {  a  [CO2])  ,  are  plotted  (curves  1  through  9  on  the 
figure) as  a  function  of  altitude.  Curve  10  is  the  sum  of  these 
contributions . 
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Figure  2-17  The  contributions  from  each  of  the  mechanisms  numbered  1 
through  9  on  table  2-4  to  the  quantity  ,  626  001-000  4.3um  earthlimb 
radiance  ,  are  plotted  (curves  1  through  9  on  the  figure) as  a  function 
of  altitude.  Curve  10  is  the  sum  of  these  contributions. 


oxygen  singlet  D  mechanism  and  what  is  actual  reality  would  show  up  in  the 
comparison  of  our  predictions  and  the  data  in  the  range  55  to  65  km  tangent 
altitude.  Our  model  for  this  mechanism  is  given  in  our  paper  (Kumer,  1977a). 

Considering  that  anything  very  much  different  from  our  model  (in  the  sense 
that  it  could  be  considerably  more  than  our  model  predicts)  would  show  up  in 

the  low  altitude  data,  and  it  does  not.  We  believe  that  our  model  is 

perfectly  adequate  for  consistent  interpretation  of  the  SPIRE  data,  and  we  see 
no  evidence  in  these  data  for  any  0(^D)  model  that  is  very  much  different 
from  ours.  We  must  add  that  there  are  a  number  of  factors  to  be  considered 
when  judging  the  applicability  of  any  model  for  the  0(1D)  mechanism.  The 
product  of  the  O(^D)  production  rate  with  the  efficiency  for  transfer  of 
excitation  of  vibration  to  N?  is  the  pertinent  quantity.  Production  of 
0(D)  depends  on  the  0^  number  density  and  in  this  case  an  inaccurate  O3 
density  (too  large  or  too  small)  could  be  compensated  for  by  an  inaccurate 
choice  (in  the  opposite  direction)  for  the  efficiency.  Our  bottom  line  here 
is  that  agreement  with  the  data  that  we  get  with  our  present  model  argues  that 

changes  to  our  present  model  for  the  0(^0)  mechanism  are  probably  not 

warranted. 

2.3.6  Vibrational  Temperatures 

Vibrational  temperatures  and  atmospheric  kinetic  temperature  for  scan  9  are 
plotted  on  Fig.  2-18.  Curve  1  is  the  atmospheric  kinetic  temperature. 

Curve  2  is  the  N2  vibrational  temperature.  Curves  J+2  designated  the 
vibrational  temperature  of  the  4.3  ym  C02  band  designated  by  J  in 
Table  2-3.  For  example  curve  3  is  the  vibrational  temperature  of  the  626 
C02  (001)  state.  Curve  8  is  the  vibrational  temperature  of  the  636  C02 
(Oil)  state,  etc.  The  vibrational  temperatures  are  defined  by 

T  =  -  1.439  v* / ln(B)  (2.3) 


where  B  is  defined  by 
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Figure  2-18  Vibrational  temperatures  and  atmospheric  kinetic 
temperature  for  scan  9  are  shown  plotted  on  the  figure  .  The  curve 
labelled  by  1  the  atmospheric  kinetic  temperature.  The  curve 
labelled  by  2  is  the  N2  vibrational  temperature.  Curves  labelled  by  J+2 
designate  the  vibrational  temperature  of  the  4.3 um  C02  band  J  in  table 
2-3  above  . 


(2.4) 


[C02  (Jvjj  =  q'B  [C02  (000)] 

and  where  v1  is  the  upper  state  energy  in  units  cm-*,  q’  is  the  degeneracy 
of  the  upper  state.  The  stanaaro  that  [xj  symbolizes  the  number  density  of 
the  species  x  is  utilized  above.  For  the  definition  of  vibrational 
temperature  is  given  by  Eq.  (2.3)  above  with  B  given  by 

[N2(V  -  1)]  -  b  [n23 

The  vibrational  temperatures  listed  on  Fig.  2-18  were  calculated  by  methods 
and  techniques  described  in  the  (1977 a)  paper  by  Kumer.  The  radiance 
predictions,  which  agree  very  well  with  the  SPIRE  data  (see  Fig.  2-11  through 
2-14  above)  as  is  discussed  above,  are  related  to  these  vibrational 
temperature  calculations  in  that  they  specify  the  excited  state  populations 
that  are  used  in  the  calculation  of  the  radiances.  We  believe  this  good 
agreement  between  our  predictions  and  the  data  provides  confirmation  of  our 
model  for  the  daytime  C02  atmospheric  emission  processes  and  confirmation 
for  our  predictions  for  daytime  C02  and  N2  vibrational  temperatures  that 
were  published  in  the  (1977a)  paper. 

2.3.7  Band  and  Line  Shape  Altitude  Dependence  in  the  Limb  Radiance  Calculation 

Referring  to  the  discussion  of  radiance  calculation  that  appears  in 
section  2.2.3  we  note  that  our  4.3  pm  model  has  not  been  completely  upgraded 
to  include  the  effect  of  altitude  dependence  in  the  line  and  band  shapes  as 
has  been  done  for  the  C02  626  010-000  band  in  our  15  urn  model.  Such  an 
upgrade  of  the  4.3  pm  model  would  involve  considerable  effort  and  was  beyond 
the  scope  of  the  effort  we  report  here.  One  of  the  problems  with  upgrading 
the  4.3  pm  model  is  that  so  many  bands  are  important,  i.e.,  the  fundamental 
band  becomes  relatively  unimportant  for  altitude  <  100  km  for  example.  We 
did,  however,  find  a  way  to  approximately  account  for  line  and  band  shape 
altitude  dependence  in  the  4.3  pm  radiance  calculation  for  all  the  bands  by 
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utilizing  an  analytic  approximation  for  the  transmission  function  matrix 
elements  that  appear  on  Fig.  2-10.  This  approximation  can  be  inferred  from 
the  work  of  Neuendorf fer  (1982)  on  developing  an  analytic  expression  for  the 
absorption  of  a  single  non-overlapping  line  for  a  given  mass  path  length  of 
absorber  gas,  and  a  mass  averaged  lover  the  path  length)  value  for  the 
pressure  broadening  parameter  and  for  the  Doppler  width  of  that  line.  The 
Neuendorffer  approximation  is  an  improvement  over  previous  analytic 
approximations  in  that  it  is  generally  applicable  through  both  the  optically 
thick  and  the  thin  regime,  and  in  that  it  may  be  utilized  for  arbitrary  ratio 
of  Lorentz  to  Doppler  broadening.  It  is  restricted  to  homogeneous  paths,  but 
can  be  applied  to  the  earthlimb  case  by  replacing  the  inhomogeneous  earthlimb 
mass  path  with  an  equivalent  homogeneous  pseudo-mass  path  specified  by  the 
inhomogeneous  earthlimb  mass  path  averaged  values  of  pressure  (Lorentz)  and 
Doppler  broadening. 

To  verify  the  Neuendorffer  approximation  and  to  illustrate  the  expected  effect 
on  the  4.3  urn  radiance  calculation  we  show  on  Fig.  2-19  several  calculations 
of  the  earthlimb  width  function  for  the  626  001-000  band  that  differ  in  the 
assumptions  that  are  employed  in  the  specification  of  band  and  line  shape. 

The  width  function  for  a  CO2  band  is  defined  on  Fig.  2-9  above.  The 
physical  significance  of  the  width  function  is  that  a  good  approximation  to 
the  earthlimb  radiance  is  given  by  the  product  of  the  width  function 
calculated  along  an  earthlimb  viewing  line  of  sight  for  a  specific  tangent 
altitude  and  the  emission  source  function  value  at  the  tangent  point.  The 
emission  source  function  is  the  quantity  s  that  is  defined  in  section  2.2. 

On  Fig.  2-19  curve  1  is  the  quantity  oN  for  earthlimb  geometry,  plotted  as  a 
function  of  tangent  altitude.  In  the  optically  thin  case  the  width  function 
is  equal  to  oN  independent  of  the  band  shape  and  of  the  value  used  for  Doppler 
and/or  Lorentz  widths;  this  is  seen  to  be  the  case  for  tangent  altitudes  > 

130  km  for  the  626  001-000  band  for  earthlimb  viewing  geometry. 

The  width  function  calculations  shown  on  Fig.  2-19  are  characterized  by 
parameter  "a,"  called  the  Lorentz  parameter,  that  is  defined  as  the  inverse  of 
the  Doppler  e-fold  1/2  width  divided  by  the  Lorentz  half  width  at  half  max. 
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The  width  function  calculations  presented  on  Fig.  2-19  differ  by  the  value  of 
"a"  that  is  used  for  each  calculation,  and  by  the  distribution  of  rotational 
levels,  which  determine  bandshape,  that  is  used  for  each  calculation.  Curve  2 
was  computed  on  the  assumption  that  a  =  10“  ^  over  the  entire  tangent 
altitude  range.  Curve  2  was  computed  by  the  methods  we  have  discussed 
previously  (Kumer  and  James,  1974).  Curves  3,  4,  and  5  were  calculated  by  the 
Neuendoeffer  approximation.  Curve  5  was  calculated  for  the  pure  Doppler  case 
i.e.,  a  =  0  over  the  entire  altitude  range;  Curve  4  was  calculated  for  the 
same  condition  as  curve  2,  i.e.,  a  =  10“^,  and  curve  3  was  computed  by 
utilizing  a  mass  averaged  altitude  dependent  value  of  "a"  that  is  plotted  on 
Fig.  2-20.  For  curves  2,  4,  and  5  the  band  and  Doppler  line  shape  are  set  to 
be  given  by  225  deg,  considered  to  be  generally  representative  of  mesopheric 
conditions.  Curve  3  is  the  most  realistic  case.  For  curve  3  the  altitude 
dependence  of  the  line  shape  is  accounted  for  by  combining  the  Neuendorffer 
analytic  approximation  with  mass  averaging.  The  altitude  band  shape 
dependence  due  to  temperature  variation  is  also  accounted  for  by  utilizing 
mass  averaging  to  assign  a  mean  temperature  to  each  limb  viewing  path  (this 
procedure  is  identical  with  that  used  to  assign  a  Doppler  width  to  each  limb 
viewing  path). 

The  agreement  between  curves  2  and  4  provides  mutual  verification  for  our 
previous  method  and  for  the  Neuendorffer  analytic  approximations.  The 
comparisons  indicate  that  there  may  not  be  a  lot  of  effect  on  the  C02  626 
001-000  4.3  um  radiance  calculation  due  to  altitude  dependence  of  line  and 
band  shape  above  approximately  80  km.  However,  below  80  km  it  would  appear 
that  there  should  be  a  large  effect  due  to  the  increasing  pressure  broadening 
with  decreasing  altitude.  One  would  expect  this  since  in  first  approximation 
the  radiance  is  given  by  the  product  of  the  width  function  and  the  emission 
source  function  S.  This  effect  is  illustrated  by  curve  3. 

A  comparison  of  the  data  for  scan  9  and  our  predictions  is  shown  on 

Fig.  2-21.  The  predictions  shown  on  Fig.  2-21  differ  from  those  shown  on 

Fig.  2-11.  For  predictions  shown  on  Fig.  2-11,  the  altitude  dependent  value 

for  "a"  as  shown  on  Fig.  2-20  was  used  in  a  Neuendorffer  like  calculation 
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Figure  2-19  The  quantity  a  nl  (curve  l)for  earthlimb  geometry  and  the 
earthlimb  viewing  width  function  W(  a  n  )  that  is  calculated  by 
utilizing  various  assumptions  for  the  line  and  band  shapes  as 
discussed  in  the  text  are  plotted  on  this  figure  as  a  function  of 
tangent  altitude. Curves  2,4  and  5  use  altitude  independent  line  and 
band  shapes, curve  3  uses  a  more  realistic  altititude  dependent 
"earthlimb  path  mass  averaged"  line  and  band  shape  . 


radiance  IN  WATTS/ C STER*CM**2 ) 


Figure  2-21  Same  as  Figure  2-11  except  in  this  case  the  calculation 
utilizes  the  same  altitude  independent  line  and  band  shapes  that  were 
utilized  above  in  the  calculation  of  the  width  function  curves  2  and  4 
on  Figure  2-20  .  By  contrast  the  the  more  realistic  "earthlimb  path 
mass  averaged”  line  and  band  shapes  were  used  for  the  calculations 
shown  on  Figure  2-11  through  Figure  2-14  . 
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accounting  for  altitude  dependent  band  and  line  shape  variation  of  the 
transmission  function  matrix  elements  (these  are  defined  on  Fig.  2-10).  This 
is  required  to  achieve  an  "exact"  limb  radiance  calculation  as  defined  on 
Fig.  2-10.  By  contrast  the  altitude  independent  band  shape,  and  line  shape 
characterized  by  a  =  10  ,  is  used  in  the  calculation  of  the  transmission 

function  matrix  elements  that  were  used  to  produce  the  predictions  that  are 
shown  on  Fig.  2-21. 

By  inspection  of  the  comparison  between  data  and  prediction  on  Fig.  2-11  with 
that  on  Fig.  2-21,  it  is  seen  that  better  agreement  can  indeed  by  obtained 
between  the  data  and  the  predictions,  all  the  way  down  to  about  53  km  tangent 
altitude,  by  use  of  the  more  realistic  calculation  which  takes  into  account 
the  altitude  dependence  of  the  band  and  line  shapes.  Below  53  km  the 
predictions  on  Fig.  2-11  rapidly  become  much  greater  than  the  data.  This 
happens  because  our  model  utilizes  the  assumptions  of  non-overlapping  lines. 
This  is  a  good  assumption  above  53  km  but  at  lower  altitudes  this  assumption 
begins  to  break  down.  For  example  the  square  equivalent  width  of  the 
strongest  line  of  the  626  001-001  band  is  about  1.3cm-*  for  the  40  km 
tangent  altitude  limb  viewing  path;  this  is  to  be  compared  with  the  average 
line  spacing  of  about  1.6  cm-*  in  that  band.  In  this  case  then  (according 
to  the  assumption  of  non-overlapping  lines)  all  of  the  lines  of  the  weaker 
bands  that  are  located  in  the  same  spectral  region  as  the  strong  line  are 
incorrectly  making  an  additional  contribution  to  the  limb  radiance  at  40  km. 

This  is  not  the  case  at  higher  altitudes  since  it  can  be  shown  that  the 
equivalent  square  width  of  the  strong  lines  varies  like  exp  (-  z/h)  where  h  is 
the  atmospheric  scale  height  (for  saturated  lines  with  strong  Lorentz  wings  as 
is  the  case  for  strong  lines  of  the  626  001-000  band  near  40  km).  In  this 
case  a  line  that  is  1.3  cm-*  wide  at  40  km  slims  out  to  about  0.13  cm-* 
equivalent  width  at  approximately  55  km.  This  is  considerably  less  than  the 
average  line  spacing  in  band.  Considering  the  number  and  relative  strength  of 
the  weaker  bands  that  are  contributing,  it  becomes  reasonable  that  the 
non-over-lapping  assumption  is  adequate  for  altitudes  >  55  km.  This  is 
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adequate  for  altitudes  755  km.  This  is  verified  to  be  the  case  by  the 
comparisons  of  our  predictions  with  the  data  that  are  shown  on  Fig.  2-11  and 
Fig.  2-21. 

On  Fig.  2-11  the  discrepancy  between  theory  and  data  below  approximately  55  km 
could  be  reduced  by  using  a  radiance  calculation  that  takes  line  overlap  into 
account.  One  way  to  do  this  might  be  to  modify  LOWTRAN  so  that  we  could 
supply  the  non-LTE  CO2  vibrational  temperature  for  the  calculation  of  the 
Planck  black  body  radiance  (which  is  required  for  calculation  of  the 
atmospheric  radiance  by  LOWTRAN)  and  at  the  same  time  maintain  the  standard 
LOWTRAN  procedure  of  utilizing  the  atmospheric  kinetic  temperature  for  the 
calculation  of  atmospheric  transmissivity  that  is  also  required  for  the 
LOWTRAN  calculation  of  the  radiance.  This  is  a  promising  approach  to 
obtaining  better  agreement  with  the  data  in  the  55  to  40  km  region  that  we 
were  not  able  to  pursue  in  the  scope  of  the  currnt  effort.  We  believe  this 
would  be  a  useful  task  to  pursue  in  future  efforts. 

2.3.8  Possibility  for  Azimuthal  SPIRE  Scan  Pattern  Misregistration 

In  section  2.3.3  we  mentioned  that  our  predictions  were  generally  slightly 
larger  than  the  data  (see  Figs.  2-11  through  2-14).  Here  we  explore  the 
possibility  that  the  solar  elevation  angles  appropriate  for  analysis  of  the 
SPIRE  data  may  be  smaller  than  is  currently  accepted  due  to  some  inadvertent 
azimuthal  misregistration  of  the  SPIRE  scan  pattern.  Since  the  directional 
data  for  SPIRE  pointing  were  somewhat  degraded,  this  possibility  does  have 
some  viability  as  pointed  out  above.  The  SPIRE  data  have  been  well  calibrated 
vertically  by  AFGL  personnel  by  using  prominent  atmospheric  infrared  features 
such  as  the  OH  maximum,  the  HNO-j  maximum,  the  C03  4.3  um  knee,  etc.  The 
azimuth  however  may  not  be  as  well  calibrated,  thus  an  azimuthal  rotation  of 
the  entire  SPIRE  scan  pattern  (the  pattern  with  currently  accepted  aximuth  is 
shown  in  page  205  of  the  AFGL  Technical  Memo.  AFGL-TM-18,  Nadile,  1979) 
towards  the  night  side  would  result  in  lower  solar  elevation  angles  for  the 
sunlit  scans  than  are  currently  accepted. 
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On  the  basis  of  the  currently  accepted  SPIRE  scan  pattern  azimuth  Nadile  and 
Grieder  predict  a  solar  elevation  angle  7.7°  for  tangent  point  of  the  earth's 
surface  for  scan  11  (the  closest  of  the  sunlit  scans  to  the  dark  region). 
However  when  we  plot  the  window  region  3.9-4. 1  pm  data  versus  tangent  altitude 
down  to  negative  tangent  altitudes  (i.e.,  for  looking  at  the  hard  earth)  it 
appears  (see  Fig.  2-22)  that  scan  11  much  more  closely  resembles  the  dark  data 
than  the  sunlit  data.  This  might  suggest  that  the  solar  elevation  at  ground 
level  for  scan  11  is  less  than  the  nominal  7.70°. 

This  same  behavior  is  not  conclusively  mirrored  in  other  window  spectral 
regions  of  the  SPIRE  data  however.  For  example,  AFGL  personnel  find  that  in 
the  2.2  ym  window  region,  the  radiance  in  the  earth  viewing  dark  scans  that  is 
below  NESR;  is  heavily  saturated  for  the  sunlit  scans  9,  10,  and  12,  but 
fairly  mildly  saturated  for  scan  11.  Since  scans  9,  10,  and  12  are  heavily 
saturated,  it  is  difficult  to  tell  whether  these  represent  full  sun  and  scan 
11  a  terminator  case  (essentially  zero  solar  elevation  angle)  or  not. 

Similarly,  AFGL  personnel  have  looked  at  the  SPIRE  data  in  the  3.65  pm  region, 
here  scan  11  is  an  order  of  magnitude  brighter  than  for  the  dark  scans.  It  is 
not  clear  that  the  3.65  pm  data  represent  reflection  of  3.65  sunlight  from  the 
earth  or  solar  pumped  non-LTE  emissions  of  CH^  and  O3  at  altitude,  (i.e., 
the  effective  earth  shadow  height  could  be  several  km  high)  so  again  it  is 
difficult  to  use  these  data  to  draw  hard  conclusions  on  what  the  solar 
elevation  angle  might  be  for  scan  11. 

Our  predictions  for  the  4.3  pm  earthlimb  radiance  for  solar  elevation  angles 
calculated  on  the  assumption  that  the  SPIRE  scan  pattern  was  in  fact  rotated 
14  deg  further  away  from  the  sun  (i.e.,  in  the  counter  clockwise  direction  as 
viewed  in  the  perspective  of  the  figure  on  page  205  of  the  Memo  AFGL-TM-18 
that  was  cited  above)  than  is  currently  believed  to  be  the  case  are  shown 
compared  against  the  data  for  scans  9  through  12  on  Fig.  2-23  through  2-26. 

It  can  be  seen  that  the  predictions  that  are  based  on  the  rotated  case  (i.e., 
those  shown  on  Figs.  2-23  through  2-26)  do  indeed  match  with  the  data  slightly 
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Figure  2-22  SPIRE  spectral  earthlimb  radiance  data  in  the  3.9  to  4.1  urn 
window  region  are  plotted  versus  tangent  altitude  h  on  this  figure. 
Note  that  scan  11  data  for  hard  earth  viewing  (ie.  for  h<0)  more 
closely  resemble  data  obtained  in  the  night  time  scans  1  through  8 
than  data  obtained  on  the  obviously  sunlit  scans  9,10  and  12. 
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Figure  2-23  Same  as  Figure  2-11  except  that  the  solar  elevation  angles 
utilized  for  the  calculation  are  reduced  by  the  assumption  that  the 
SPIRE  scan  pattern  has  been  rotated  14  degrees  further  away  from  the 
sun  than  is  currently  accepted. 


RADIANCE  IN  WATTS/ t STER*CH**2  ) 


5 


.  TANGENT  AL  T  I  TUOE  IN  KM 


Figure  2-24  Same  as  Figure  2-12  except  that  the  solar  elevation  angles 
utilized  for  the  calculation  are  reduced  by  the  assumption  that  the 
SPIRE  scan  pattern  has  been  rotated  14  degrees  further  away  from  the 
sun  than  is  currently  accepted. 
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as  Figure  2-13  except  that  the  solar  elevation  angle 

sPTRP^oran^n  H?e  CaiCUl?tion  are  reduced  by  the  assumption  that  the 
SpiRE  scan  pattern  has  been  rotated  14  degrees  further  away  from  the 
sun  than  is  currently  accepted. 
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Figure  2-26  Same  as  Figure  2-14  except  that  the  solar  elevation  angles 
utilized  for  the  calculation  are  reduced  by  the  assumption  that  the 
SPIRE  scan  pattern  has  been  rotated  14  degrees  further  away  from  the 
sun  than  is  currently  accepted. 


better  than  do  the  predictions  generated  on  the  basis  of  no  rotation  from  the 
accepted  scan  azimuth  (shown  in  Figs.  2-11  through  2-14  above). 

Considering  both  the  "nadir  4.0  pm  window"  data  and  the  high  altitude  daytime 

4.3  pm  data,  we  do  seem  to  have  reason  to  suggest  that  there  may  be  some 
residual  uncertainty  in  the  registration  of  the  azimuth  of  the  SPIRE  scan 
patterns.  We  are  reluctant  to  push  this  suggestion  very  hard  though,  since  in 
neither  case  is  the  argument  conclusive  as  much  as  it  is  suggestive.  That  is, 
in  either  case  alternate  possibilities  that  might  explain  the  evidence  that  we 
have  construed  in  these  cases  to  suggest  azimuthal  misregistration  can  be 
invoked.  For  example,  in  section  2.3.3  we  cite  several  mechanisms  alternate 
to  azimuthal  misregistration  which  could  contribute  to  the  overprediction  that 
is  generally  evident  in  Figs.  2-11  through  2-14.  The  latter  predictions  are 
based  on  the  currently  accepted  azimuthal  registration.  Similarly,  we  have 
discussed  in  this  section  some  of  the  effects  that  render  the  "nadir  4.0  pm 
window"  data  more  suggestive  than  conclusive,  that  there  could  be  some 
misregistration  in  the  SPIRE  scan  pattern  azimuth. 

2.4  ANALYSIS  OF  THE  NIGHT  TIME  SPIRE  4.3  pm  DATA 
2.4.1  Introduction 

In  this  section  we  discuss  a  brief  correlative  analysis  of  the  2.7  and  4.3  pm 
data  obtained  in  the  eight  night  time  SPIRE  data  scans  and  a  more  detailed 
quantitative  analysis  of  a  representative  sample  of  the  night  time  SPIRE 
4.3  pm  data,  namely  quantitative  analysis  of  the  night  time  scans  2  and  8. 

The  scope  of  our  effort  limited  us  to  quantitative  analysis  of  just  a  portion 
of  the  data.  Scans  8  and  2  are  the  extreme  "south"  and  "north"  scans 
respectively,  thus  we  thought  these  scans  might  provide  extremes  of  night  time 
variability.  Scan  1  was  slightly  further  "north"  than  scan  2  but  we  did  not 
choose' to  analyze  this  scan  in  detail  since  it  was  not  complete,  i.e.,  it  did 
not  extend  up  to  high  enough  altitude.  The  data  were  processed  as  described 
in  section  2.3.1  in  an  attempt  to  extract  effects  such  as  drift  of  the  DC 
level  in  the  data. 


The  comparison  of  our  predictions  for  the  night  time  4.3  ym  earthlimb  against 
the  data  for  scans  2  and  8  generally  supports  the  position  in  our  paper  (Kumer 
et  al.,  1978)  that  night  time  CO,,  4.3  pm  emission  in  the  altitude  region 
>  -  75  km  is  dominated  by  collisional  transfer  of  vibrational  excitation  from 
0H(v  >  0)  to  N2>  followed  by  subsequent  transfer  to  CO2  and  multiple 
scattering  in  the  4.3  pm  bands.  The  correlative  study  of  all  eight  scans  did 
not  conclusively  support  or  refute  this  position. 

We  also  find  that  additional  mechanisms  are  required  to  achieve  a  better  fit 
to  the  data.  The  requirement  for  these  additional  mechanisms  becomes 
particularly  noticeable  above  approximately  90  to  95  km.  We  looked 
quantitatively  into  the  possibility  of  two  mechanisms  to  explain  the 
additional  high  altitude  excitation,  (1)  the  direct  excitation  by  sunlight, 
and  (2)  an  excitation  related  to  the  recombination  of  oxygen,  i.e.,  a 
Chapman-like  mechanism.  Neither  of  these  can  account  for  the  extra  excitation 
mechanism  that  is  required  at  high  altitude.  The  major  problem  with  the 
Chapman  mechanism  is  that  it  requires  an  efficiency  for  producing  N2 
vibrational  excitation  from  the  energy  released  by  the  exothermic  reaction 
0  +  0  +  M  -  >  O^+M  that  is  nearly  200  percent.  We  discuss  results  obtained 
and  requirements  on  this  Chapman-like  mechanism  for  the  extra  source  of  C02 
(v3)  excitation  above  90  to  ~  95  km  in  detail  in  section  2.4.2. 

We  must  caution  that  these  findings  should  be  considered  as  preliminary  since 
we  did  not  have  the  resources  in  the  scope  of  the  effort  that  we  are  reporting 
to  produce  a  scientifically  complete  investigation.  To  do  this  would  require 
us  to  quantitatively  analyze  all  eight  data  scans  that  were  obtained  in  night 
time  conditions.  Since  the  two  mechanisms  cited  above  could  not  explain  the 
extra  excitation  above  90  to  95  km,  it  is  required  that  we  should  consider  a 
number  of  additional  mechanisms.  Some  candidate  mechanisms  that  need  to  be 
considered,  but  were  not  in  the  scope  of  the  effort  we  are  reporting  on,  would 
include  precipitation,  transport  of  C02  (vj)  by  multiple  scattering  of 
sunlight  from  the  adjacent  sunlit  region  of  the  atmosphere,  or  a  broader 
source  of  0H(v)  than  we  actually  used,  one  that  extends  to  higher  altitudes 
than  currently  envisaged  for  production  of  0H(v). 
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2.4.2  Comparison  of  Predictions  and  Data  for  Scans  2  and  8 

On  Figs.  2-27  through  2-29  we  compare  the  data  for  scan  2  with  our 
predictions.  As  above,  the  data  are  represented  by  curve  1,  the  predictions 
are  curve  2,  the  626  CC^  001-000  band  contribution  to  the  predictions  is 
curve  3,  and  the  contribution  of  the  weak  bands  is  curve  4.  We  progressively 
add  in  mechanisms  to  achieve  a  better  fit  in  going  from  Figs.  2-27  to  2-29. 

The  mechanisms  we  consider  for  the  night  time  case  are  given  in  Table  2-5. 

Note  that  absorption  of  sunlight  is  not  included  in  Table  2-5.  Since  the  hard 
earth  shadow  heights  were  less  than  140  km  for  many  of  the  "dark"  scans  we 
specifically  investigated  the  posibility  that  absorption  of  sunlight  could 
play  a  role  in  producing  the  extra  source  of  excitation  of  CO.,  (V3)  that 
is  required  at  high  altitude.  We  found  that  this  was  not  an  important  factor 
since  the  large  absorption  along  the  long  slant  paths  through  the  atmosphere 
renders  the  effective  shadow  height  for  absorption  of  sunlight  by  CO^  to  be 
approximately  80  to  110  km  higher  than  the  hard  earth  shadow  height. 

The  non-airglow  associated  mechanisms,  thermal  collisions  and  absorption  of 
4.3  um  earthshine,  are  utilized  to  attempt  to  model  the  aata  on  Fig.  2-27. 
These  are  mechanisms  1,  4,  and  5  listed  in  Table  2-5.  It  can  be  seen  that 
these  mechanisms  do  not  produce  sufficient  excitation  to  adequately  model  the 
data. 

On  Fig.  2-28  we  add  in  the  OH  airglow  related  mechanism  2  in  Table  2-5.  This 
mechanism  has  been  discussed  by  Kumer  et  al  (1978)  where  it  was  shown  that 
excitation  transfer  via  mechanism  2  in  Table  2-5  of  about  0.12  ergs/(S-cm2) 
was  required  to  achieve  agreement  with  a  set  of  quiet  night  time  4.3  um  zenith 
radiance  data  obtained  4-11-74  by  AFGL  and  Utah  State  University  personnel. 

To  achieve  the  modeling  agreement  shown  by  Fig.  2-28  required  transfer  from 
0H(V  >  0)  via  mechanism  2  in  Table  2-5  of  0.19  ergs/(S-cm2)  in  the  half 
width  altitude  range  83  to  89  km.  It  can  be  seen  that  an  additional  source  of 
excitation  above  90  km  is  required  to  achieve  a  better  fit  to  the  data. 
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Figure  2-27  The  4.3 urn  data  for  dark  scan  2  are  given  by  curve  1.  The 
predictions  are  labelled  by  curve  2 /components  due  to  the  CO2  626 
001-000  band  by  curve  3  and  due  to  the  "weak"  bands  by  curve  4.  The 
non-airglow  associated  mechanisms  1,4  and  5  in  Table  Figure  2-5  are 
used  to  caculate  the  predictions. 


The  predictions  we  show  on  Fig.  2-29  include  a  second  postulated  airglow 
related  mechanism.  The  "Chapman-like"  mechanism  is  shown  by  item  3  in 
Table  2-5.  The  excitation  Ex  in  units  {vibrational  quanta)/(S  -  cm  )  due 
to  this  postulated  mechanism  is  given  as  a  function  of  altitude  by 

Ex  =  Ef  [M]  [0]2 

and  the  value  Ef  =  1.1  E  -  31  cm^/S  is  required  to  achieve  the  quality  of 
fit  shown  on  Fig.  2-29.  Since  the  rate  constant  k  =  3  x  E  -  33  cm  /S  for 
the  reaction 

0+0+M->02+M  ( 

and  since  there  are  5ev  released  by  reaction  (2.5),  we  are  requiring  that 
about  200  percent  of  the  available  energy  finds  its  way  into  vibrational 
excitation  of  N2.  This  appears  to  exceed  the  realm  of  credibility. 


Table  2-5  MECHANISMS  FOR  NIGHT  TIME  EXCITATION  OF  C02  v3  IN  THE 
UPPER  ATMOSPHERE 


Item 

Mechanism 

1 

Thermal  collisions  with  atmosphere  species  other  than  N2* 

2 

03+H->0H  (V>0)  +0? 

OH  (V>0)  +N2->0H  (V-l)  +  N2  (V=l) 
n2  (V=l)  +co2  <=>  N2  +  C02  (v3) 

C02  (v3)  <=>  C02  +  hv  (4.3pm) 

3 

0+0+M->02+M  +energy 

4 

Thermal  collisions  with  N2 

5 

Absorpiton  of  4.3  pm  earthshine  by  atmospheric  C02 

*  The  mechanisms  1,  4  and  5  above  have  been  discussed  in  detail  by  Kumer  and 
James  (1974),  the  mechanism  2  has  been  discussed  in  detail  in  the  paper  by 
Kumer  et  al.  (1978),  and  the  mechanism  3  is  discussed  in  detail  in  the 
present  report. 


On  Fig.  2-30  we  show  the  comparison  of  our  predictions  with  scan  8  data.  Both 
the  OH  and  the  Chapman-like  mechanism  were  utilized  to  provide  the  fit  as 
shown.  Again  the  transfer  of  0.19  ergs/(S-cm^)  from  OH  (V  >  0)  was  used  to 
achieve  the  fit  to  the  scan  8  data  that  is  shown  in  Fig.  2-30.  In  order  to 
invoke  the  Chapman-like  mechanism  to  achieve  the  quality  of  fit  as  shown 
required  =  -  3.4E-31  cm®/S.  This  differs  from  the  value  required  to 
fit  the  scan  2  data  by  a  factor  -  3.  This  provides  further  evidence  that  the 
Chapman-like  mechanism  is  not  a  good  candidate  for  explaining  the  extra  energy 
source  for  two  reasons;  (1)  the  requirement  on  is  inconsistent  from  scan 
2  to  scan  8,  and  (2)  the  value  required  for  Ef  for  scan  8  is  so  large  that 
it  becomes  inconveivable  that  the  Chapman-like  mechanism  can  supply  the 
required  energy. 

2.4.3  Night  Time  Vibrational  Temperatures 

The  atmospheric  kinetic  temperature  and  the  vibrational  temperatures  of  N2 
(V  =  1)  and  C02  (001)  which  are  calculated  for  scan  2  on  the  basis  of  the 
non-airglow  associated  mechanisms  1,  4,  and  5  in  Table  2-5  are  shown  plotted 
as  a  function  of  altitude  on  Fig.  2-31.  Changes  to  the  vibrational 
temperatures  obtained  by  progressively  adding  in  the  mechanisms  2  and  3  in 
Table  2-5  are  illustrated  by  Fig.  2-32  and  2-33. 

2.4.4  Comparison  of  the  SPIRE  2.7  and  4.3  urn  Night  Time  Earthlimb  Data 

In  order  to  determine  if  we  could  establish  a  correlation  between  the  2.7  and 
4.3-pm  night  time  earthlimb  on  a  scan  by  scan  basis,  to  support  the  hypothesis 
that  the  OH  airglow  mechanism  is  the  dominant  excitation  mechanism  for  both 
these  spectral  regions  in  the  night  time  earthlimb,  we  performed  a  brief 
examination  of  the  2.7  and  4.3  pm  earthlimb  in  the  altitude  region  of  the  OH 
airglow  emission  layer  in  the  2.5  to  3.1  pm  and  in  the  4.1  to  4.5  pm  bands  for 
all  eight  "dark"  scans.  The  spectral  radiances  in  the  2.7  pm  region  varied 
between  1.3  and  2.3  x  10”^  watts/  (sq.  cm-sr-um)  and  the  spectral  radiances 
in  the  4.3  pm  region  varied  from  6  to  10  x  10~7  watts/(sq.  cm-sr-um). 
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Figure  2-31  Atmospheric  kinetic  temperature (curve  1) ,N2  and  CC2(001) 
vibrational  temperatures  (curves  2  and  3  respectively)  are  shown 
here .Mechanisms  are  those  utilized  to  generate  the  predictions  for 
scan  2  that  are  shown  on  Figure  2-27. 
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Figure  2-32  Atmospheric  kinetic  temperature (curve  1),N2  and  CO?(001) 
vibrational  temperatures  (curves  2  and  3  respectively)  are  shown 
here .Mechanisms  are  those  utilized  to  generate  the  predictions  for 
scan  2  that  are  shown  on  Figure  2-28. 


73 


f,2"3?  AtroosPheric  kinetic  temperature (curve  1) ,N2  and  C02 (001) 

K^faii0I!aliteraperature8  <curves  2  and  3  respectively)  are  shown 

o^»!!*'>e£v!a!?iSmS  a5e  th0Se  utilized  to  generate  the  predictions  for 
scan  2  that  are  shown  on  Figure  2-29. 


First,  there  appeared  to  be  only  a  small  correlation,  if  any  at  all,  on  going 
from  one  scan  to  the  next  in  either  of  the  wavelength  regions.  The  noise  and 
the  DC  drift  in  the  baseline  of  the  spectral  data  could  have  contributed  to 
difficulty  in  establishing  any  readily  apparent  "one  scan  to  the  next 
correlation"  in  a  given  band.  Also  the  altitude  resolution  is  approximately 
equal  to  the  half  width  of  the  OH  emission  layer  and  this  could  also  have 
introduced  some  problems  in  trying  to  identify  any  apparent  trend  in  the  data 
on  going  from  "one  scan  to  the  next". 

Secondly,  we  did  not  find  any  readily  identifiable  correlation  between  the  2.7 
and  4.3  pm  data  on  a  scan  by  scan  comparison  basis.  However,  correlation  on  a 
scan  by  scan  basis  in  the  2.7  and  4.3  pm  data  is  difficult  to  identify  for  the 
same  reasons  given  above  for  the  difficulties  faced  in  identifying  any 
correlation  or  trend  within  a  given  band  on  going  from  one  scan  to  the  next. 
Also,  the  multiple  scattering  of  photons  by  C02  would  tend  to  smooth  spatial 
structure  that  may  be  present  in  the  OH  airglow  source.  Since  the  relaxation 
time  constant  for  C02  emission  is  6  to  8  minutes  in  this  altitude  region, 
this  too  would  tend  to  smooth  temporal  structure  in  the  OH  airglow  source. 

We  conclude  that  the  comparison  study  between  the  2.7  and  the  4.3  pm  night 
time  earthlimb  did  not  show  any  readily  apparent  correlation  on  a  scan  by  scan 
basis,  but  for  reasons  given  above  this  should  not  be  taken  to  mean  that  the 
OH  mechanism  is  not  the  dominant  mechanism  for  producing  C02  (v^)  in  the 
altitude  region  80  to  90  km. 

2.5  ANALYSIS  OF  DAYTIME  C02  CONTRIBUTION  TO  THE  2.7  pm  EARTHLIMB 

2.5.1  Estimated  C02  Contribution  to  the  Sunlit  2.7  pm  Earthlimb 

In  our  January  1981  presentation  at  the  AIAA  meeting  held  in  St.  Louis  (Kumer, 
1981)  we  showed  that  absorption  of  sunlight  at  2.7  pm  by  the  626  C02  species 
followed  by  reemission  at  2.7  pm  could  not  produce  sufficient  radiance  to 
explain  the  2.7  pm  daytime  SPIRE  data.  This  stimulated  work  by  Sharma  to 
investigate  the  possibility  that  absorption  of  sunlight  by  C02  at  ~2.0  and 
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1.6  pm  followed  by  fluorescence  at  2.7  pm  could  be  an  important  contributor  to 
the  SPIRE  sunlight  2.7  pm  earthlimb.  Likewise  Degges  investigated  the 
prospect  that  resonance  scattering  of  sunlight  at  2.7  pm  by  might  also 

be  an  important  contributor. 

In  the  event  Sharma  and  Degges  showed  that  both  of  these  mechanisms  were 
indeed  very  important  in  the  formation  of  the  2.7pm  earthlimb.  Additionally 
we  looked  into  the  prospect  that  absorption  and  reemission  of  2.7  pm  sunlight 
by  the  minor  isotopic  C02  species  636,  628,  and  627  could  also  play  a 
significant  role  in  the  formation  of  the  2.7  pm  sunlit  earthlimb,  and  found 
indeed  that  it  did.  By  considering  the  contributions  of  all  these  mechanisms 
Sharma  and  Nadile  (1981)  showed  that  the  2.7  pm  sunlit  SPIRE  data  could  be 
reasonably  well  understood. 

On  Fig.  2-34  we  show  our  estimated  prediction  for  the  contribution  to  the 

2.7  pm  earthlimb  that  is  due  to  absorption  and  reemission  of  2.7  pm  sunlight 
by  C02.  Curve  2  is  the  component  that  is  due  to  absorption  and  reemission 
of  2.7  pm  sunlight  by  the  major  isotopic  species  626  C02.  The  sum  of  all 
these  components  is  curve  1  on  Fig.  2-34.  The  contribution  from  absorption 
and  subsequent  reemission  of  2.7  pm  sunlight  by  the  636,628  and  627  isotopic 
species  of  C02  is  given  by  curve  3. 

The  predictions  that  we  show  for  curves  2  and  3  are  a  factor  2  smaller  than  we 
had  reported  previously  in  our  March  through  April  1981  progress  report. 

Since  the  time  of  that  work  we  have  found  and  corrected  a  very  troublesome  bug 
in  our  program  that  had  introduced  an  overprediction  by  a  factor  of  2  in  our 
curves  2  and  3  on  Fig.  2-34.  We  found  and  corrected  this  error  just  recently 
in  the  process  of  analyzing  the  daytime  4.3  pm  earthlimb  data  as  is  described 
above  in  section  2.3.  Our  program  had  originally  been  functioning  accurately 
but  we  introduced  the  bug  when  we  made  an  attempt  to  deal  with  the  effect  of 
the  reaction 
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Figure  2-34  Predictions  for  the  CO2  contribution  (curve  1)  to  the 
sunlit  2.7  wm  earthlimb  that  results  from  the  absorption  of  2.1m 
sunlight  by  CO2  .followed  by  subsequent  re-emission  .are  shown  here. 
Curve  2  is  the  component  due  to  the  626  CO2  specie  and  curve  3  is  the 
component  due  to  the  minor  isotopic  species  of  CO2  . 
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We  were  concerned  with  this  reaction  because  the  states  with  =  2  are 
forbidden  to  radiate  in  the  2.7  pm  band,  thus  the  effect  of  this  reaction 
would  be  to  reduce  the  predictions  cn  Fig.  2-34  by  another  factor  2  if  the 
rate  constant  for  reaction  2.6  is  as  large  as,  or  comparable  to,  the  rate 
constant  for  the  reaction 

C02(v3)  +  N2  -  >  C02  +  N2  { V=l)  (2.7) 

Since  our  calculations  on  Fig.  2-34  represent  an  underprediction,  one  might 
infer  that  the  rate  constant  for  the  reaction  2.6  is  not  as  large  as  that  for 
the  reaction  2.7.  We  had  hoped  to  measure  the  rate  constant  for  reaction  2.6 
in  our  complimentary  laboratory  experiment  but  circumstances  prevented  this. 
Details  are  provided  in  section  3.0  below  which  deals  with  our  complimentary 
laboratory  effort.  Future  effort  should  be  directed  at  development  of  an 
accurate  method  to  account  for  reaction  2.6  quantitatively  in  the  calculation 
of  its  effects  on  the  2.7  pm  earthlimb,  ana  to  establish  realistic  estimates 
of  the  magnitude  of  the  rate  constant,  i.e.,  it  is  not  of  much  interest  unless 
it  is  comparable  to  or  larger  than  for  reaction  2.7. 

2.5.2  Preliminary  Efforts  to  Spectrally  Model  the  2.7  pm  Earthlimb 

The  2.7  pm  sunlit  SPIRE  data  lend  themselves  to  spectral  analysis  which  may 
possibly  yield  fairly  precisely  the  contribution  from  each  of  the  candidate 
sources.  To  this  end  we  began  a  preliminary  effort  to  adapt  a  very  general 
and  comprehensive  software  package  that  has  been  developed  partially  under  the 
Lockheed  Independent  Research  Program  and  partially  under  the  ongoing  contract 
NAS5-26173  for  the  purpose  of  calculation  of  spectral  limb  radiance  and  for 
the  decomposition  of  spectral  earthlimb  data  into  the  components  that  result 
from  various  species  and/or  mechanisms.  The  package  includes  capability  for 
relatively  rapid  calculation  at  spectral  resolution  av  >  vp  where  vq  is 
the  Doppler  width  of  the  radiating  molecular  species.  Arbitrary  slit 
functions  may  be  utilized.  The  calculations  may  be  performed  for  any  of 
20  molecular  species  compiled  on  the  latest  version  of  the  AFGL 
McClatchey/Rothman  tape.  On  Fig.  2-35  we  show  a  sample  calculation  of 
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Figure  2-35  The  contribution  to  the  sunlit  2.7um  spectral  earthlimb 
that  is  due  to  absorption  and  re-emission  of  sunlight  by  CO2  is  shown 
here.  The  simulated  data  are  calculated  for  resolution  and  sensitivity 
that  approximates  SPIRE  performance. 


earthlimb  spectral  radiance  at  70  km  tangent  altitude  due  to  absorption  of 
sunlight  at  2.7  pm  by  CO2  followed  by  reemission  at  2.7  pm.  The  calculation 
was  performed  for  conditions  of  solar  illumination  approximately  applicable  to 
sunlit  scan  9  of  the  SPIRE  experiment. 

A  facility  of  this  software  is  that  the  contribution  to  spectral  radiance  due 
to  any  bands  of  interest  to  the  programmer  may  be  calculated  and  displayed 
independently.  The  software  also  has  facility  to  use  these  components  as  base 
functions  to  be  fit  to  actual  data,  thereby  providing  the  quantitative  means 
to  decompose  the  SPIRE  data  into  the  components  that  result  from  various 
species  and/or  mechanisms.  The  software  may  also  be  employed  in  an  inversion 
mode  to  determine  either  the  number  density  of  the  emitting  species,  or,  in 
some  cases,  to  determine  the  vibrational  temperature  of  the  upper  levels  of 
the  transition.  The  projected  capability  of  this  improved  software  package 
which  is  being  developed  under  the  programs  cited  above  differs  from  our 
previous  capability  in  that  (1)  spectral  earthlimb  radiance  is  calculated, 
rather  than  molecular  band  integrated  earthlimb  radiance,  (2)  the  capability 
may  be  applied  to  all  20  species  on  the  latest  Rothman  update  of  the 
McClatchey  tape,  (3)  the  new  capability  does  not  utilize  physical  models  to 
calculate  vibrational  temperature;  instead,  the  vibrational  temperatures  must 
be  (a)  input  from  an  independent  calculation  such  as  can  be  provided  in 
certain  cases  by  our  previous  software  capability,  or  (b)  derived  by  an 
inversion  process  directly  from  earthlimb  data  such  as  the  SPIRE  data,  or  data 
that  are  to  be  obtained  in  the  ELIAS,  CIRRIS,  and  SPIRIT  experiments. 

A  logical  development  for  the  improved  software  package  would  be  the  detailed 
adaptation  to  the  analysis  of  DNA/AF6L  earthlimb  problems  that  require  this 
type  of  detailed  spectral  earthlimb  modeling  capability.  This  would  include 
further  detailed  application  against  the  SPIRE  data,  to  glean  the  added 
information  that  is  contained  in  the  spectral  signature  of  the  SPIRE  data. 

Such  an  effort  need  not  be  limited  to  the  2.7  pm  region  but  in  fact  would  be 
useful  throughout  the  1.4  to  16.5  pm  spectral  region.  The  improved  software 
may  also  be  used  for  analysis  of  data  to  be  obtained  in  upcoming  earthlimb 
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experiments  such  as  SPIRIT  and  CIRRIS  that  are  tc  employ  very  high  spectral 
resolution. 


2.6  HORIZONTAL  STRUCTURE  CONSIDERATIONS 


To  this  point  the  discussion  in  section  2  has  been  concentrated  on  the 
analysis  of  the  SPIRE  data.  In  this  section  we  digress  from  this  theme,  to 
utilize  the  SPIRE  data  at  15  iim  to  assist  in  constructing  a  reasonably 
realistic  scenario  for  the  occurrence  of  horizontal  structure  in  earthlimb 
radiance  for  "long  period"  gravity  waves  (approximately  30  min  or  more). 


There  are  cases  in  the  literature  (Noxon,  1978;  or  Faire  and  Murphy,  1972  are 
a  small  cross  section  of  these)  which  suggest  that  gravity  waves  with 
atmospheric  temperature  perturbations  of  the  order  40°K  are  not  uncommon  in 
the  mesopause  region.  Also,  there  is  evidence  that  suggests  that  horizontal 
wavelengths  may  be  comparable  with  length  scales  of  interest  to  the 
community.  For  example,  Peterson  (1979)  discusses  a  bright  structured  visual 
OH  airglow  event  suggestive  of  gravity  waves  with  16  km  horizontal  wavelengths 
which  occured  on  4  and  5  April  1978.  On  the  other  hand,  gravity  waves  ducted 
through  the  mesopause  region  may  have  horizontal  wavelengths  of  hundreds  of  km. 


A  plausible  approximation  for  a  gravity  wave  perturbation  to  the  model 
atmosphere  (shown  on  Fig.  2-2  and  2-3  above)  that  we  used  for  modeling  the 
SPIRE  data,  as  is  discussed  in  section  2.2,  is  given  by  equations  (2.8)  and 
(2.9)  below. 


T  =  T0(Z)  jl  +  A(z)  sin  [k2(z  -  80  km)  +  kx  X0  +  4/~]  J 

ni  =  n0i(7)  jl  -  A (z)  sin  [kz(z  -  80  km)  +  kx  X0  +  #r]  j 

Here  Tq  an(j  nQi  are  the  temperature  and  number  densities  of  the 
unperturbed  atmosphere. 


(2.8) 

(2.9) 


For  this  sample  case  we  use  kx  =  2ir/140  km,  kz  =  2u/20  km,  w  =  3.06  x 
10~3  s'1,  =  0.08  (km/sec)3,  g  =  0.0095  km/s3,  H  =  6  km  and  y  = 

7/5  and  =  (y  -  1 ) g^/C2  =  4.5  x  10'4  s'3  so  that  (from  Hines 
Equations  (33)  through  (36))  we  get 

«)  =  «g(kx/k2)  =  (2.12  x  10'2)  (20/140)  =  3.03  x  10'3  s'1 

Thus  the  period  is  35  min,  considerably  slower  than  the  Brunt  frequency.  For 
this  case  it  can  be  shown  that  the  approximate  formulation  of  Hines  (1960)  can 
be  used  to  show  that  the  relative  phase  and  amplitudes  of  the  fractional 
perturbations  in  temperature  and  number  density  are  as  shown  in  Eqs.  2.8  and 
2.9  above. 

The  fractional  amplitude  A(z)  of  the  perturbation  is  selected  so  that  the 
temperature  perturbation  is  ~  40  deg  at  -  95  km.  The  specific  altitude 
dependence  of  the  fractional  perturbation  is  given  by: 

A  =  0.056  exp  [(z  -  80)/2H] 

for  z  <_  100  km  where  z  is  altitude  in  km  and  H  =  6  km.  But  for  z  >  100  km  we 
introduce  an  empirical  damping  fator,  so  in  this  case 

A  =  0.056  exp  [(z  -  80)/2H]  exp  [(100  -  z)/H]. 

To  illustrate  the  possible  infrared  horizontal  structure  effect  of  a  plausible 
approximation  for  a  gravity  wave  perturbation  to  the  model  atmosphere  (shown 
on  Fig.  2-2  and  2-3)  that  we  used  for  modeling  the  SPIRE  data  as  is  discussed 
in  section  2.2,  we  show  earthlimb  radiance  profiles  due  to  C0£  010  »  000 
emission  in  the  14  to  16  uni  band  on  Fig.  2-36.  Curve  D  is  representative  of 
the  SPIRE  data  as  shown  on  Tig.  2-1.  Curve  U  is  calculated  via  our  most 
recent  techniques  as  described  above  in  section  2.2.  Curve  U  is  identical 
therefore  to  our  predictions  for  the  15  ym  earthlimb  that  are  compared  with 
the  SPIRE  15  ym  data  on  Fig.  2-6.  It  is  calculated  for  the  case  of  the 
unperturbed  model  atmosphere  we  have  used  for  analysis  of  the  SPIRE  data  that 


« 


82 


2-36  Barthlimb  radiance  in  the  14  to  16  um  region  for  various  cases  as 
explained  in  tne  text  are  shown  here  .Curve  D  is  represenative  of  the 
SPIRE  data .Curve  0  is  a  calculation  based  on  the  unperturbed  atmospheric 
model (figures  2-2  and  2-3  above)  that  we  have  utilized  in  our  SPIRE  data 
analysis  effort.  Curves  0  and  n  are  calculations  based  on  gravity  wave 
perturbations  of  the  SPIRE  data  analysis  model  atmosphere  that  differ  by 
n  in  phase . 
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is  illustrated  by  Fig.  2-2  and  2-3.  Curves  0  and  *  apply  for  calculations 
that  are  performed  with  the  atmosphere  perturbed  by  a  gravity  wave  as  given  by 
Eqs.  (2.8)  and  (2.9)  above  for  the  phase  angle  \lr=  0  and  u  respectively. 

Note  that  this  example  assumes  plane  gravity  waves  that  are  propagating  in  a 
direction  that  is  perpendicular  to  an  earth  limb  sensor  line  of  sight  (LOS), 
hence  there  is  no  gravity  wave  induced  horizontal  variation  along  the  line  of 
sight  for  this  example.  The  gravity  wave  model  that  we  have  utilized  is  a 
simple  variation  of  the  isothermal  Hines  (1960)  model  and  is  useful  for  this 
illustration.  However,  more  realisitic  calculations  of  the  gravity  wave 
perturbation  of  temperature  and  density  in  a  real  atmosphere  are  necessary  to 
assess  the  effect  of  gravity  waves  more  reliably  than  we  have  done  here.  This 
would  require  utilization  of  the  latest  results  of  Tuan  and  coworkers  (Private 
Communication,  Tuan,  1982). 

2.7  IMPROVEMENTS  SUGGESTED  IN  THE  CURRENT  HIGH  ALTITUDE  INFRARED  RADIANCE  MODEL 

The  SPIRE  experiment  and  subsequent  analysis  provides  a  vehicle  to  verify  and 
upgrade  the  AFGL  High  Altitude  Radiance  Model.  This  is  an  extremely  useful 
thing  to  do  since  experience  has  shown  that  derivatives  of  the  High  Altitude 
Radiance  Model,  such  as  the  OSC  Code  earthlimb  model,  are  widely  utilized  in 
industry  for  evaluating  sensor  performance  in  sensor  simulation  studies.  Too 
frequently  the  codes  are  applied  by  personnel  with  little  knowledge  of  the 
applicability  of  the  code  to  their  particular  problem,  and  not  much  of  an  idea 
as  to  when  the  code  may  not  be  applicable  to  their  particular  problem.  Since 
this  kind  of  utilization  is  sometimes  practiced  (certainly  not  the  recommended 
practice  by  any  means)  it  is  necessary  that  the  possible  shortcomings  of  the 
model  be  explicitly  spelled  out  so  that  misuse  of  the  model  is  minimized,  and 
that  the  model  be  kept  up  to  date  to  the  extent  that  the  availability  of  data, 
such  as  the  SPIRE  data,  and  the  analysis  thereof  permits. 
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In  this  discussion  we  address  our  remarks  to  the  version  of  the  Current  High 
Altitude  Radiance  Model  (CHARM)  that  is  described  by  T.C.  Degges  in  the  report 
of  Humphrey  et  al.  (1981).  Major  ingredients  of  the  CHARM  include  the 
(1)  physical  radiative  transfer  model,  (2)  non-LTE  emission  mechanisms, 

(3)  the  identification  of  the  spectral  bands  of  a  given  radiating  specie  that 
must  be  included  in  CHARM  in  order  to  perform  adequate  calculations,  and 

(4)  the  temperature  and  composition  of  the  input  model  atmosphere.  Our 
discussion  will  address  items  (1)  through  (3). 

Assumptions  that  go  into  the  physical  radiative  transfer  model  that  are  used 
in  CHARM  are  listed  on  page  12  of  the  report  by  Humphrey  et  al.  (1981).  The 
analysis  of  the  14  to  16-ym  radiance  that  we  report  in  section  2.2  above 
indicates  that  utilization  of  the  assumption  (6)  "temperature  variations  at 
different  altitudes  do  not  affect  the  radiative  transfer"  in  the  CHARM  does 
indeed  introduce  serious  error  in  the  altitude  region  above  95  km.  In 
section  2.2  we  have  shown  that  the  inclusion  of  band  and  line  shape  variation 
versus  altitude,  which  results  directly  from  "temperature  variations  at 
different  altitudes,"  is  one  of  the  major  effects  which  must  be  accounted  for 
to  improve  on  pre-SPIRE  modeling  of  the  15-ym  SPIRE  data  is  achieved. 

We  also  believe  that  the  CHARM  could  be  improved  by  utilizing  line  strengths 
that  are  properly  weighted  across  the  band  as  we  described  above  in 
section  2.2,  rather  than  by  use  of  the  current  CHARM  assumption  (4)  "line 
strengths  are  used  as  if  each  line  were  at  the  band  center,  neglecting 
wavelength  variation  across  the  band."  The  method  we  use  to  account  for 
variation  of  line  strength  across  the  band  has  been  discussed  by  Kumer  and 
James  (1974)  and  is  also  discussed  in  section  2.2  above.  We  believe  this 
method  could  be  included  into  CHARM  and  would  result  in  an  improved  version  of 
CHARM  that  would  suffer  little  or  no  increased  complexity  and  run  time. 

The  analysis  of  SPIRE  data  we  report  here  has  also  given  some  insight  as  to 
additional  non-LTE  mechanisms  that  are  required  to  improve  the  CHARM. 

Inclusion  of  the  faster  vtr  rate  constant  for  quenching  of  C02  (010)  by 
atomic  oxygen  at  mesospheric  temperatures  is  required  to  better  model  the 


15— ym  earthlimb.  This  was  the  second  effect  that  we  found  in  section  2.2 
above  to  be  very  important  to  achieve  an  adequate  model  of  the  15-ym  SPIRE 
data.  The  inclusion  of  the  0(*D)  mechanism  is  required  to  achieve  a  better 
model  for  the  daytime  4.3-ym  earthlimb.  For  the  nighttime  4.3-ym  earthlimb 
our  analysis  of  the  SPIRE  data  in  section  2.4  above  has  shown  that  the 
transfer  of  energy  from  0H(V  >  0)  is  required  to  be  included  into  CHARM  in 
order  to  adequately  model  the  data  in  the  approximate  altitude  range  75  to 
90  km.  The  analysis  also  indicates  that  another,  as  yet  unidentified 
mechanism  needs  to  be  specified  in  order  to  account  for  the  underprediction  of 
SPIRE  nighttime  4.3-ym  data  above  approximately  90  km  that  we  discussed  in 
section  2.4  above. 

The  version  of  CHARM  that  is  described  in  the  report  by  Humphrey  et  al.  (1981) 
does  not  include  the  "weak"  bands  2  through  8  that  we  have  presented  in 
Table  2-3  above.  However,  we  have  shown  in  section  2.3  above  that  these  bands 
are  major  contributors  to  both  the  day  and  nighttime  4.3-ym  earthlimb  in  the 
altitude  region  from  approximately  60  to  105  km  and  are  required  to  achieve  an 
accurate  model  of  the  4.3-ym  SPIRE  data.  Hence  an  improved  and  realistic 
version  of  the  CHARM  should  include  these  bands.  The  analysis  also  shows  that 
a  larger  number  of  bands  than  presently  included  in  CHARM  should  in  fact  be 
included  in  order  to  adequately  model  the  15-ym  earthlimb.  Likewise  Sharma 
has  shown  by  analysis  of  the  SPIRE  2.7-ym  data  that  the  bands  which  account 
for  the  fluorescence  of  sunlight  at  2.0  and  1.6  ym  should  be  included  in  the 
CHARM  in  order  to  adequately  model  the  day  time  2.7-ym  SPIRE  data.  Finally, 
bands  of  the  minor  isotopic  species  of  CO2  have  been  shown  to  be  required  in 
all  three  regions  of  the  SPIRE  data  that  we  have  subjected  to  detailed 
analysis,  namely  15,  4.3,  and  2.7  ym. 

In  this  section  we  have  summarized  a  number  of  aspects  of  the  CHARM  that  the 
SPIRE  data,  and  the  analysis  thereof,  show  require  improvement.  By  the  time 
our  report  is  published  we  would  expect  that  some  of  these  deficiencies  would 

have  been  corrected  in  the  CHARM.  j 
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Section  3 

LABORATORY  STUDY  OF  C02  FLUORESCENCE 
3.1  EXPERIMENTAL  OBJECTIVES  AND  APPROACH 

During  the  time  period  covered  in  this  report  we  undertook  a  number  of 


experiments  aimed  at  studying  the  processes: 

C02  (000)  +  hv  (2.7  ym)— C02  (101) 

(3-la) 

co2  (101)  +  co2  (000)— co2  (100)  +  co2  (001) 

(3-lb) 

and 

C02  (001)— C02(000)  +  hv  (4.3  ym) 

( 3— 1C ) 

C02  (101)— C02(000)  +  hv  (2.7  ym) 

( 3— Id ) 

C02  (101)— C02(100)  +  hv  (4.3  ym) 

(3-le) 

The  primary  goal  of  these  experiments  was  to  obtain  experimental  values  for 
the  branching  ratio  which  described  what  fraction  of  C09  in  the  101  level 
decayed  via  process  (3-le)  (emission  at  4.3  ym)  as  compared  with  process 
( 3— Id )  (emission  at  2.7  ym). 

The  approach  was  to  produce  molecules  in  the  101  level  via  process  (3-la)  by 
using  a  tunable  2.7  ym  laser  to  excite  the  C02.  By  comparison  of  the 
relative  intensity  of  emission  at  4.3  ym  and  2.7  ym  the  branching  ratio  can  be 
determined. 

There  are  a  number  of  factors  which  complicate  these  measurements  and  which 
make  the  use  of  a  pulsed  laser  a  desirable  technique  for  obtaining  the  desired 
results. 

In  the  absence  of  any  emission  measurements  it  would  be  expected  that  the 
branching  ratio  could  be  calculated  from  known  values  of  the  absorption  band 
strengths  for  the  4.3  ym  band  and  the  2.7  ym  band.  While  this  is  true,  there 
a  number  of  reasons  for  testing  an  emission  measurement  aside  from  the  value 
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of  having  a  direct  confirmation  of  the  branching  ratio  to  compare  with  the 
calculated  ratio.  The  branching  ratio  is  affected  by  many  collisional 
processes  which  result  in  a  transfer  of  the  initial  vibrational  excitation  to 
other  molecules  or  else  cause  a  redistribution  of  that  energy  among  various 
vibrational  modes. 

Depending  on  the  number  of  densities  of  the  colliding  partners,  many  of  these 
processes  occur  on  time  scales  smaller  than  the  radiative  lifetime.  As  a 
result,  the  actual  branching  ratios  obtained  under  a  given  set  of  conditions 
may  differ  considerably  from  the  values  expected  for  the  isolated  molecules  in 
a  hypothetical  collision  free  environment. 

The  single  most  important  process  which  affects  the  laboratory  measurements 
and  leads  to  considerable  difficulty  in  the  measurement  is  process  (3-lb)  (the 
sharing  of  the  excitation  between  two  molecules  as  a  result  of  collision  of 
excited  C02  with  C02  in  the  ground  state).  While  for  the  isolated 
molecules  the  lifetime  for  radiating  at  2.7  pm  is  longer  than  that 
corresponding  to  4.3  urn,  it  turns  out  that  the  exchange  shown  in  (3-lb)  occurs 
with  such  a  high  efficiency  that  any  2.7  pm  radiation  is  quenched  at  a  rate 
fast  compared  to  emission  at  4.3  pm.  This  requires  that  the  2.7  pm 
measurement  be  made  on  a  time  scale  which  is  much  smaller  than  that  required 
at  4.3  pm.  The  emission  at  4.3  pm  is  not  affected  by  this  because  emission  at 
4.3  pm  from  the  states  (101)  and  (001)  presumably  takes  place  at  the  same  rate 
and  therefore  process  (3-lb)  which  redistributes  the  excitational  energy 
should  not  affect  the  overall  rate  of  emission  at  4.3  pm. 

The  rate  constant  for  reaction  (3-lb)  is  1.3  x  10“^°  cm^  s_*.  (Finzi 
and  Moore,  1975).  Therefore  at  a  pressure  of  the  order  of  10~2  torr,  for 
example,  the  lifetime  of  the  101  level  is  of  the  order  of  21  pS,  which  is  more 
than  two  orders  of  magnitude  smaller  than  the  4.3  pm  band  radiative  lifetime. 
This  clearly  points  to  the  need  to  be  able  to  measure  the  relative  emission  on 
a  time  scale  much  faster  than  required  for  pure  radiative  decay. 
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In  actual  measurements,  a  pressure  of  10-2  torr  would  be  unsuitable  because 
under  the  conditions  of  our  experiment  which  utililzed  a  small  fluorescence 
cell,  the  results  would  be  effected  by  deactivation  as  a  result  of  wall 
collisions.  In  our  experiments  we  utilized  C02  in  the  presence  of  50  torr 
of  Argon.  Under  these  conditions  the  Vj  vibration  responsible  for 
4.3  emission,  is  quenched  with  a  lifetime  of  approximately  400  ys.  Still,  the 
2.7  ym  radiation  is  quenched  by  reaction  3-16  at  a  much  more  rapid  rate  even 

_p 

at  C02  pressures  as  low  as  10  torr. 

The  need  for  using  a  laser  in  order  to  study  the  branching  ratio  becomes 
apparent  from  the  above  considerations  of  the  erergy  transfer  process  (3-lb). 
In  order  to  have  a  reasonable  rate  for  the  transfer  reaciton,  small  number 
densities  are  required.  At  the  small  number  densities,  the  use  of  a 
continuous  source  would  not  provide  enough  energy  for  easy  measurement  of  the 
resultant  signal.  Thus  a  laser  source  which  can  provide  large  initial 
excitation  of  the  (101)  level  is  required.  Furthermore,  a  pulsed  laser  is 
ideal  because  of  the  rapid  decay  of  the  (101)  level  following  initial 
excitation. 

Previous  preliminary  experiments  which  are  described  in  the  final  report 
"Fluorescence  Experiment  and  Auroral  Data  Evaluation  to  Improve  Prediction  of 
Nuclear  Atmosphere  Infrared  Background"  DNA  4906F  by  Kumer  and  James  had 
indicated  that  our  chromatix  laser  which  was  tunable  at  2.7  ym  would  be 
adequate  for  this  purpose. 

In  the  following  section  we  describe  our  experimental  effort  and  some  of  the 
modifications  to  our  original  experiments  which  were  made  in  an  attempt  to 
overcome  certain  experimental  difficulties.  Prior  to  this  discussion, 
however,  we  first  give  a  brief  summary  of  the  types  of  measurements  which  we 
expect  to  be  able  to  measure  with  a  sufficiently  stable  pulsed  laser  source. 

Equations  (3-la)  through  (3-le)  represent  only  a  portion  of  the  processes 
which  are  of  importance  in  atmospheric  applications.  In  addition  to  the  band 
(101),  there  is  another  band  which  is  excited  by  radiation  near  2.7  ym.  That 
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is  the  (02°1)  band  at  2.7  pm.  Therefore  all  of  the  reactions  shown  in 
equations  (3-1)  will  also  occur  with  the  (101)  level  replaced  by  the  (02°1) 
band. 

The  levels  (02°1)  and  (101)  of  C02  are  quite  closely  spaced  (about  100 
cm--*-)  and  as  a  result  can  become  equilibrated  via  collisions  of  the  type: 

C02  (10*1)  +  M — ►  C02  (02*1)  +  M  (3-2) 

Therefore  any  inital  excitation  of  the  101  level  should  also  result  in  the 
(02°1)  level  becoming  populated.  The  rate  constant  for  this  reaction  can  be 
estimated  as  about  3  x  10  crrrs  (Sharma,  1968)  which  is  the  value 
obtained  for  the  corresponding  reaction  with  v3  =  0. 

In  our  laboratory  measurements  with  50  torr  of  Argon  present  in  the  cell,  this 
equilibration  should  occur  on  a  time  scale  much  faster  than  reaction  (3-lb)  so 
that  regardless  of  whether  the  (10“l)  band  or  the  ( 02 ° 1 )  band  is  excited,  the 
emission  observed  near  2.7pm  will  include  the  contributions  of  both  bands: 

( 10°  1 )  - •'(000)  +  hv2_7  (3-3) 

(02*1) - -(000)  +  hv2  7? 

Also  fluorescence  near  4.3  um  will  include  bands 

(101) - ^  (10*0)  +  hv4>3  (3-4) 

(02*1) - *-(02°0)  +  hv4’3 

The  reactions  (3.2)  would  affect  the  branching  ratio  measurements  in  that  the 
value  obtained  would  represent  an  average  value  for  the  2.7  pm  and  2.77  pm 
bands. 

In  addition  to  the  states  (10°1)  and  (02°1)  which  are  split  by  Fermi 
Resonance,  there  is  another  nearby  state  (0221)  which  differs  from  (02°1)  in 
having  vibrational  angular  momentum.  This  t  state  with  i  =  2  would  also  be 
expected  to  be  formed  in  collisions  of  the  type. 
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(3-5) 


co2  ((J^}}  +  M-C02  (0221)  +  M 

The  state  (02^1)  can  still  radiate  at  4.3  pm  to  produce  (0220);  however, 
the  transition 

(0221) — -hv  +  (000)  (3-6) 

is  forbidden  with  the  result  that  process  (3-5)  would  tend  to  reduce  the 
intensity  in  the  2.7  pm  region  and  also  increase  the  apparent  lifetime.  This 
process  would  therefore  affect  the  measured  value  of  the  branching  ratio.  The 
exact  amount  by  which  the  branching  ratio  would  be  affected  would  depend  on 
whether  or  not  process  (3-5),  which  involves  a  change  a  Z=  2,  has  the  same 
rate  constant  as  process  (3-2).  Later  in  this  report  we  suggest  ways  in  which 
this  might  be  determined. 

A  major  process  affecting  excited  CO^  in  the  atmosphere  and  one  which  can 
easily  be  studied  in  the  laboratory  is  the  transfer  of  excitation  between 
CO^  and  N2  according  to: 

C02  (V1,V2,V3  +  1)  +  N2X-C02(V1’V2’V3>  +  N2  (v=1)  (3-7) 

This  process  of  rapid  sharing  of  vibrational  excitation  between  C02  and  N2 
affects  the  time  dependence  of  the  C02  emission  in  the  atmosphere.  In  the 
laboratory,  information  concerning  this  ratio  can  be  obtained  from  the  study 
of  the  time  dependence  of  the  C02  emission  as  different  amounts  of  N2  are 
mixed  with  the  C02  in  the  sample  cell.  Preliminary  measurements  involving 
this  process  are  discussed  later  in  this  report. 

3.2  SUMMARY  OF  EXPERIMENTAL  MEASUREMENTS 

In  our  report  (DNA  4906F)  we  described  preliminary  measurements  of  the  4. 3/2. 7 
branching  ratio  in  which  a  boxcar  integrator  was  used  to  measure  the  2.7  pm 
signal.  One  of  the  major  difficulties  encountered  in  these  earlier 
measurements  was  the  fact  that  there  was  a  large  scattered  pulse  from  the 
2.7  pm  laser  light  used  to  excite  C02.  In  those  earlier  experiments  the 


preamplifier  used  to  measure  the  scattered  2.7  pm  light  had  a  response  time  of 
the  order  of  10  microseconds.  This  resulted  in  a  large  wide  (>10ps)  pulse 
from  the  exciting  laser  light  pulse  which  overlapped  and  masked  the  2.7  pm 
reemission  from  CO2.  When  this  signal  was  amplified  it  saturated  the 
amplifier.  To  avoid  this  probelem  we  worked  at  low  CO?  pressures  (5  x 
10  torr)  to  increase  the  lifetime  of  the  101  level  and  placed  the  boxcar 
window  a  long  way  away  (70  ps)  from  the  exciting  pulse  in  order  to  eliminate 
the  interference  between  this  exciting  line  and  the  resonance  fluorescence  at 
2.7  pm.  The  low  pressures  and  the  fact  that  we  were  many  half-lives  away  from 
the  peak  of  the  fluorescence  required  sacrificing  signal  intensity. 

In  an  attempt  to  remedy  this  situation  we  built  an  amplifier  with  a  response 
of  1.5  ps,  expecting  that  we  could  work  at  higher  pressures  with  a 
corresponding  increase  in  signal,  and  also  separate  the  scattered  pulse  from 
the  fluorescent  pulse. 

As  it  turned  out,  we  continued  to  have  problems  with  measuring  the  signal  at 
2.7  pm. 


Our  earlier  experiments  utilized  a  cylindrical  fluorescence  cell  in  which  the 
two  detectors  (one  for  4.3  pm  and  one  for  2.7  pm  fluorescence)  viewed  along 
the  axis  of  the  cylinder  through  two  CaF^  windows.  A  2.7  pm  filter  was 
placed  in  front  of  one  detector  and  a  4.3  pm  filrer  was  placed  in  front  of  the 
other.  The  2.7  pm  laser  pulse  was  directed  perpendicular  to  the  axis  through 
two  windows  which  had  been  placed  in  holes  drilled  through  the  cylinder 
walls.  In  the  earlier  experiments  the  laser  light  made  just  one  pass  through 
the  cell.  The  geometry  of  this  is  shown  in  Fig.  3-1.  In  order  to  enhance  the 
signals  we  built  a  multiple  pass  cell  by  placing  parallel  mirrors  inside  the 
cell  to  give  the  geometry  shown  in  Fig.  3-2.  This  cell  design  was  chosen 
because  it  was  one  which  would  increase  the  number  of  passes  without  producing 
a  significant  loss  in  solid  angle  between  the  excitation  region  and  the 
detectors.  In  most  of  our  subsequent  experiments  the  mirrors  were  positioned 
to  give  seven  passes  through  the  cell. 


The  large  amplified  scattered  pulse  originating  from  the  exciting  laser  light 
continued  to  create  problems  with  our  electronics.  This  was  due  in  part  to 
saturation  of  electronics  and  equally  as  much  due  to  the  fact  that  our  laser 
was  extremely  unstable.  For  each  exciting  pulse  a  corresponding  fluorescence 
pulse  as  monitored  by  4.3  um  fluorescence  did  not  always  occur.  The  4.3  um 
pulses  which  did  occur  were  stored  in  a  Tracor-Northern  signal  averager  so 
that  after  a  run  of  several  minutes  (with  5  laser  pulses  per  second)  the 
4.3  um  signal  was  well  defined.  Since  the  signal  at  2.7  um  always  contained 
the  effects  of  the  tail  of  the  2.7  um  scattered  laser  exciting  line  plus 
(presumably)  some  occasional  contribution  from  resonance  fluorescence  at 
2.7  um  results  obtained  with  the  boxcar  integrator  were  unstable. 

During  the  course  of  several  months  we  attempted  a  number  of  further 
modifications  to  our  system  in  an  attempt  to  overcome  these  difficulties. 

None  of  these  were  satisfactory  in  leading  reproducible  results  for  the  ratio 
of  signal  strengths  obtained  in  the  2.7  um  and  4.3  um  channels. 

A  brief  summary  of  the  types  of  modifications  that  we  made  is  given  in  the 
following  section. 

Our  first  attempt  to  eliminate  the  problem  arising  from  the  scattering  of  the 
initial  2.7  um  laser  pulse  was  to  construct  a  differential  amplifier  and  look 
directly  at  the  initial  exciting  pulse  and  compare  this  with  the  2.7  um  pulse 
of  readiation  scattered  in  the  fluorescense  cell.  By  adjusting  the  signal 
from  those  two  sources  we  expected  to  have  zero  output  from  the  differential 
amplifier  when  no  fluorescence  was  taking  place.  When  fluorescence  did  occur, 
we  expected  that  the  longer  lived  (longer  than  the  initial  excitation  pulse) 
2.7  um  fluorescence  signal  would  appear  at  the  output.  Again  this  was  not 
completely  satisfactory  because  of  saturation  and  corresponding  long  recovery 
times,  also  it  was  difficult  to  accurately  match  the  signal  from  the  two 
outputs  of  either  differential  amplifier. 
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Fig.  3-1  Experimental  Arrangement  for  Measurement  of  Fluorescence  at  4.3pm 
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Fig.  3-2  Multiple  Pass  Cell 


Continuing  efforts  involved  rejecting  the  differential  amplifier  and  using  the 
actual  4.3  ym  signal  itself  to  trigger  the  signal  averager  and  boxcar 
integrator  which  recorded  the  4.3  ym  fluorescence  and  2.7  ym  re-radiation.  By 
using  a  voltage  comparator  we  could  trigger  the  signal  averager  and  boxcar 
only  when  the  fluorescence  at  4.3  ym  was  greater  than  a  preset  value.  Since 
we  were  not  observing  fluorescence  with  every  pulse  of  the  laser,  this 
eliminated  averaging  scattered  pulses  at  2.7  ym  which  did  not  include  a 
fluorescence  component.  It  also  served  to  improve  signal  to  noise  ratio  in 
the  4.3  ym  channel  since  pulses  less  than  some  pre-determined  value  were  not 
included.  Still  we  continued  to  have  irreproducible  results  in  the  2.7  ym 
channel . 

Several  further  modifications  included  the  use  of  fast  switching  whereby  the 
signal  at  2.7  ym  was  not  admitted  to  the  amplifying  electronics  until  several 
microseconds  after  the  initial  pulse  so  as  to  eliminate  saturation  of 
electronics  caused  by  this  source.  When  these  modifications  failed  to  yield 
reproducible  results  we  finally  replaced  the  boxcar  integrator  which  appeared 
to  be  unstable  with  an  integrating  circuit  built  specifically  to  cover  the 
time  intervals  of  interest  for  measuring  fluorescence  at  2.7  ym.  A  diagram  of 
this  configuration  is  shown  in  Fig.  3-3. 

In  spite  of  the  numerous  modifications  to  our  system  we  were  not  able  to 
obtain  reproducible  results  for  the  branching  ratio  measurements.  During  the 
time  period  when  these  various  modifications  were  being  made  we  observed  that 
the  chromatix  laser  was  quite  erratic  in  its  output  and  that  this  effect 
seemed  to  become  worse  as  time  went  by. 

In  fact,  we  found  on  measuring  the  infrared  outout  from  the  heated  LiNbO^ 
crystal  that  wavelength  shifts  of  100  or  so  wavenumbers  would  occur  without 
any  change  in  the  control  setting.  In  retrospect  we  realize  that  most  of  our 
earlier  problems  associated  with  trying  to  measure  the  2.7  ym  fluorescence 
were  due  to  the  fact  that  our  laser  was  extremely  unstable.  When  we  realized 
that  our  laser  was  not  putting  out  a  constant  frequency  we  returned  the  system 
to  chromatix  for  repairs.  As  it  turned  out  a  transistor  in  the  heater  circuit 
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Fig.  3-3  Geometry  and  Electronics  for  Two-Channel  Run 


which  maintained  the  LiNbOj  crystal  at  a  constant  temperature  had  failed 
which  ultimately  resulted  in  fluctuations  of  many  wavenumbers  from  day  to 
day.  In  our  earlier  report,  DNA  4906F,  we  had  mentioned  that  the  2.7  y m 
filter  placed  in  front  of  the  detector  rejected  much  of  the  exciting  laser 
light  but  transmitted  a  significant  portion  of  the  101  and  021  bands.  Thus 
originally  the  interference  of  scattered  laser  light  was  not  as  damaging  as  it 
was  in  our  later  experiments,  since  in  these  later  experiments  the  exciting 
2.7  urn  was  changing  in  frequency  from  day  to  day  by  an  amount  sufficient  to 
cause  it  to  fall  within  the  bandpass  of  the  2.7  yin  filter  at  times. 

The  manufacturer's  sales  personnel  had  claimed  that  their  dye  laser  plus 
associated  optical  parametric  oscillator  (heated  LiNbO-^  crystal)  would 
produce  a  wavelength  stability  of  the  order  of  0.001  cnT*  over  several 
hours.  Our  experience  was  that  this  was  not  the  case  based  on  the  fact  that 
there  were  pulse  to  pulse  variations  in  the  magnitude  of  the  fluorescence 
observed  at  4.3  ym  and  also  because  a  considerable  fraction  (sometimes  as  much 
as  70  percent)  of  the  exciting  pulses  did  not  result  in  any  fluorescence  at 
all. 

In  spite  of  the  fact  that  we  were  unable  to  obtain  reproducible  signals  at 
2.7  ym  we  could  collect  many  4.3  ym  fluorescent  pulses  and  obtain  a  decay 
curve  for  4.3  ym  emission  with  good  signal  to  noise.  This  was  due  to  the  fact 
that  the  4.3  ym  signal  is  more  intense  than  the  2.7  ym  signal  and  also  because 
the  4.3  ym  filter  rejected  scattering  of  the  2.7  ym  exciting  pulse  so  that 
there  was  not  the  problem  of  interference  and  saturation  of  electronics  that 
we  had  experienced  in  the  2.7  ym  channel. 

3.3  RESULTS  FOR  VV  BETWEEN  N2(V=1)  AND  13C1602. 

In  our  earlier  report  we  examined  the  rate  of  reaction  (3-7)  involving  the 
transfer  of  energy  between  C02  and  N2.  We  obtained  results  which  were  in 
agreement  with  the  literature  value  of  -  6  x  10-*3  cm3sec-*  for  the  rate 
constant.  In  our  present  set  of  experiments  we  were  interested  in  comparing 
the  rate  of  decay  of  4.3  ym  emission  in  the  presence  of  N2  for  both 
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C02  and  C02.  Kyy  measurement  for  reaction  (3-7)  is  important 
for  some  of  the  minor  isotopic  species  such  as  13C02  for  accurate 
predictions  of  the  C02  4.3  pm  relaxation  time  in  the  x-ray  deposition  region 
of  the  atmosphere  near  80-km  altitude.  Kumer  (1977b)  has  shown  that  the  value 
of  this  relaxation  time  can  vary  between  7  and  20  minutes  near  80  km  depending 
on  the  values  of  kvv  for  various  weak  isotopic  C02  4.3  pm  bands. 

We  made  several  measurements  of  the  quenching  of  *3C02  and  *3C02  by 
argon  and  nitrogen.  The  purpose  of  these  measurements  was  to  measure  the 
difference  in  the  transfer  rate  for  the  exchange  of  vibrational  energy  between 
C02  and  N2  for  the  different  isotopic  species  of  C02.  There  was  some 
problem  with  the  experimental  data  in  that  the  quenching  rates  appeared  to  be 
much  faster  than  our  expectations  for  both  cases.  In  the  case  of  a 
measurement  of  the  quenching  of  C02  by  argon  with  no  nitrogen  present,  we 
obtained  a  value  greater  than  the  known  rate.  We  obtained  a  lifetime  for  the 
excited  C02  of  156  pS  whereas  the  expected  value  under  our  experimental 
conditions  should  have  been  about  400  pS  (Yardley  and  Moore  1968).  This 
suggests  that  there  was  probably  some  source  of  contamination  in  the 
laboratory  sample  cell  which  was  resulting  in  a  faster  quenching  rate  and  a 
corresponding  decrease  in  the  lifetime  of  C02. 

In  spite  of  this  apparent  discrepancy  we  have  attempted  to  make  some  estimates 

as  to  the  relative  quenching  rates  of  ^C02  and  *3C02  by  N^.  To 

make  this  comparison  we  have  assumed  that  in  all  of  our  quenching  runs,  which 

used  the  same  amount  of  Argon,  that  the  same  quenching  rate  of  156  ps  would 

have  been  obtained  in  the  absence  of  N2.  what  this  amounts  to  is  the 

assumption  that  the  impurities  which  must  have  been  present  were  constant  for 

all  of  the  runs.  Based  on  this  assumption,  we  have  then  comapred  the 
12  13 

C02  and  C02  results.  To  make  this  comparison  we  have  calculated 
curves  of  the  emission  rate  as  a  function  of  time  for  various  values  of  the 
N2  transfer  rate. 


Given  that  in  the  presence  of  argon  for  our  experimental  conditions,  the 
excited  CO2*  appears  to  have  a  lifetime  of  156  ys,  we  have: 

C02*  +  (A  +  Impurities)  —  C02  +  (A  +  Impurities)  (3-8) 

-d(C02*)/dt  =  (C02*)  K  (3-9) 

where 

K  =  1/  r  =  6410  s'1 

In  addition  to  this  process,  the  excited  C02*  participates  in  the  energy 
exchange  process: 


(3-10) 


From  this  process  we  have: 

-d(C02*)/dt  .  kf(C02*)(N2)  -  kr(C02)(N2*) 


and 


-d(N2*)/dt  =  Kr(C02)(N2*)  -  ky(C02*)(N2) 

In  our  experiments  we  had  (C02)  =  (N2).  If  we  make  a  simplifying 

assumption  that  kf  =  kr  because  of  the  near  resonant  transfer,  then  we 
have: 

-d(C02*)/dt  =  k(C02*)  -  k(N2*)  +  K  (C02*) 

-d(N2*)/dt  =  k(N2* )  -  k(C02*)  (3-12) 
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where  k  is  now  given  by: 


k  a  kf(N2)  a  kr(C02) 


(3-13) 


In  the  case  of  C09»  the  measured  value  of  k*  is  reported  to  be 

1*3  3  1  C  T 

6  x  10  cm  s  (Taylor  and  Bitterman  1964).  Therefore  since  for  our 
experiments  (C02)  =  0.3  torr,  we  expect  a  transfer  rate  of  k  =  6,371  s'1. 


The  solution  of  Eqs.  (3-9)  and  (3-12)  subject  to  the  condition  that  at  time 
t  *  0  that  (N2*)  *  0  is: 


(co2)  -  (co2)t=0  ~ — 


e'Kt/2  r  f 

TT  L(Jx'K) 


-(M* 


(3-14) 


where 


vTx  = 


In  Fig.  (3-4)  and  (3-5)  we  show  calculated  curves  of  the  rate  of  emission  of 
4.3  pm  radiation  from  excited  C02*  as  a  function  of  time.  These  curves  are 
calculated  for  a  value  of  K  given  by  Eq.  (3-9)  and  various  values  of  k  of 
Eq.  (3-13).  The  experimental  results  are  then  superimposed  on  these  curves  in 
order  to  determine  the  value  of  k  which  seems  to  give  the  best  fit  of  the 
data.  From  Fig.  (3-4)  it  is  seen  that  the  data  appears  to  agree  best  with  a 
value  of  k  =  6,371  s'1  in  agreement  with  our  expectation  for  12C02.  In 
the  case  of  Fig.  (3-5),  some  of  the  data  appears  to  lie  closer  to  the  curve 
with  k  *  3186  s'1  whereas  one  of  the  runs  was  a  better  match  with  the 
k  -  1,593  s'1. 


From  the  comparison  given  above,  it  appears  likely  that  the  rate  of  exchange 

of  vibrational  energy  between  C0?  and  N?  slower  by  a  factor  of  2  or  3  for 

13  12  L  c 

C02  as  compared  with  C02.  Unfortunately  we  cannot  place  a  high 

degree  of  certainty  on  these  results  due  to  the  fact  that  we  know  from  our 
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Fig.  3 


TIME  fos) 


4  Plot  of  the  Rate  of  Emission  of  4.3  pm  Radiation  as  a  Function  of 
Time.  The  Solid  Lines  are  Calculated  from  Eq.  (7).  The  Circles 
Represent  Experimental  Data  from  Two  Runs. 
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100  200  300  400 


TIME  fos) 


Fig.  3-5  Plot  of  the  emission  of  4.3  ym  radiation  as  a  function  of  time. 

The  solid  lines  are  calculated  from  Eq.  (7).  The  Circles  and 
Triangles  Represent  Experimental  Oata  From  Different  Runs. 
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results  with  pure  Argon  that  there  must  be  some  impurities  in  the  system. 

Also  this  comparison  does  not  adequately  analyze  the  departure  of  the 
experimental  results  from  the  calculated  curves  at  small  times.  This  is 
presumeably  related  to  the  overall  temporal  response  of  the  electronic 
recording  systems. 

3.4  CONCLUSIONS  ANO  SUGGESTIONS  FOR  FUTURE  STUDIES 

Our  laboratory  studies  of  the  Fluorescence  of  CO2  at  4.3  pm  when  excited  by 
a  2.7  pm  laser  pulse  have  pointed  to  a  number  of  useful  studies  which  can  be 
made  using  this  technique.  Our  experience  also  points  to  a  number  of  problems 
which  can  arise  in  studies  of  this  type.  In  our  attempts  at  measuring  the 
branching  ratio  for  decay  of  levels  such  as  (101)  via  4.3  pm  emission  and 
re-emission  at  2.7  pm  we  encountered  many  problems  associated  with  the  need  to 
separate  the  initial  scattered  2.7  pm  pulse  from  2.7  pm  resonance  fluorescence 
following  this  excitation  pulse.  By  suitable  choice  of  infrared  filters  to 
isolate  the  2.7  pm  laser  pulse  from  most  of  the  reemitted  fluorescence  and 
fast  electronics  that  would  rapidly  recover  from  a  large  burst  of  scattered 
2.7  pm  radiation,  this  separation  should  be  achievable.  Our  efforts  to  do 
this  were  unsatisfactory  because  of  the  instability  of  our  laser  at  2.7  pm. 

The  initial  exciting  pulse  at  2.7  pm  did  not  always  occur  at  the  same 
wavelength  which  made  it  difficult  to  separate  the  2.7  pm  scattered  pulse  from 
subsequent  resonance  fluorescence  at  2.7  pm.  Any  attempts  in  the  future  to 
study  the  branching  ration  would  have  to  be  aware  of  this  problem  and  we  would 
stress  the  need  for  a  stable  2.7  pm  laser  pulse  which  could  be  tuned  to  a 
wavelength  in  the  2.7  pm  bands  which  could  be  isolated  from  much  of  the 
subsequent  band  emission  via  sharp  cutoff  infrared  filters. 

The  study  of  4.3  pm  emission  produced  via  reactions  such  as  (3-la)  and  (3-lb) 
is  much  simpler  due  to  the  fact  that  interference  from  the  2.7  pm  laser  pulse 
is  not  a  problem.  From  a  study  of  the  time  dependence  of  the  decay  of 
emission  at  4.3  p  information  can  be  gained  concerning  the  quenching  of  the 
(001)  level  by  various  species  and  also  the  rate  of  transfer  of  energy  between 
excited  CO2  and  other  species  such  as  N2.  In  order  to  obtain  accurate 
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results  for  rate  constants  from  studies  of  this  type  it  is  necessary  to  have 
precise  capability  of  measuring  pressures  and  mixing  ratios  of  gases  in  the 
pressure  ranges  from  several  microns  to  several  torr.  Contamination  by  small 
amounts  of  unwanted  gases  can  completely  destroy  a  set  of  measurements. 

Earlier  in  this  report  we  have  mentioned  that  the  initial  excitation  of  the 
levels  (101)  and  (02*1)  results  in  population  of  the  levels  (100)  and  (02°0) 
as  a  result  of  4.3  um  emission  and  also  as  a  result  of  transfer  of  vj 
excitation  between  excited  and  ground  state  C0?.  In  addition  to  these 

p 

levels  the  level  (02  0)  is  also  expected  to  be  formed  because  of  collisions 
such  as  shown  in  Eq.  (3-5)  followed  by  4.3  pm  emission  and  also  by  similar 
collisions  involving  (100)  and  (02*0)  with  various  collision  partners.  These 
processes  should  result  in  non-LTE  populations  of  levels  (100),  (02*0),  and 
(0220)  with  resultant  fluorescent  pulses  at  13.8  pm,  16.2  pm,  and  14.98  pm 
as  these  states  radiatively  decay  to  the  (01*0)  level. 

A  study  of  the  relative  intensity  and  time  dependence  of  the  emission  from 
these  three  hands  should  provide  information  on  the  rate  at  which  energy  is 
redistributed  between  these  various  levels.  The  band  at  14.98  pm  is  closely 
intermixed  with  the  (01*0)  (000)  band  at  15  mm  so  that  it  will  be  more 
difficult  to  study  than  the  13.8  and  16.2  pm  bands.  It  may  prove  feasible  to 
isolate  selected  lines  in  the  14.98  pm  band  by  a  suitable  etalon. 

Alternatively,  a  filter  consisting  of  C0?  in  an  absorption  cell  may  prove 

^  pi 

useful  in  distinguishing  between  radiation  from  the  (02  0)  (01  0)  band 

and  the  (01*0)  (000)  band. 

Based  on  our  experience  with  the  Chromatix  laser  which  did  not  perform 
according  to  the  manufacturer's  specifications  we  would  caution  those 
interested  in  pursuing  investigations  of  this  type  to  make  certain  that  they 
have  a  stable  laser  source  which  has  a  constant  unchanging  wavelength. 


A  pulsed  laser  at  2.7  pm  does  provide  an  excellent  method  for  studying  many 
fluorescent  processes  and  related  collisional  quenching  and  energy  transfer 
processes  in  CO^  which  are  of  importance  in  relation  to  the  atmospheric 
radiance  backgrounds  occurring  at  many  infrared  wavelengths. 
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2  Cy  ATTN 

LKD,  R.  Narcisi 

2  cy  ATTN 

OPR,  R.  Nadille 

DEPARTMENT  OF  THE  ARMY 

2  Cy  ATTN:  OPR,  R.  Murphey 

4  cy  ATTN:  CA,  A.  Stair 

US  Army  Electronics  RAD  Cmd 

3  cy  ATTN:  DELAS-EO,  F.  Niles 

Air  Force  Ofc  of  Scientific  Rsch 

ATTN 

AFOSR/NC 

BMD  Advanced  Technology  Ctr 

ATTN:  ATC-O.  Vi.  Davies 

Headquarters 

ATTN:  ATC-T,  M.  Capps 

Air  Force  Systems  Cmd 

ATTN 

SDR 

BMD  Systems  Cmd 

ATTN 

DLTW 

2  cy  ATTN:  BMDSC-HH 

ATTN 

DLXP 

ATTN 

OLS 

Deputy  Chief  of  Staff  for  Rsch  Dev  A  Acq 

ATTN 

DLAE 

ATTN:  DAMA-CSS-N 

Air  Force  Tech  Applications  Ctr 

US  Army  Ballistic  Rsch  Labs 

ATTN 

TD 

ATTN:  DRDAR-BLB,  J.  Hester 

ATTN 

TF 

ATTN:  DRDAR-BLA-S 

ATTN 

Tech  Library 

US  Army  Combat  Surv  A  Target  Acq  Lab 

Air  Force  Weapons  Lab 

ATTN:  OELCS-K,  S.  Kronenberg 

ATTN 

SUL 

US  Army  Foreign  Science  A  Tech  Ctr 

Air  University  Library 

ATTN:  DRXST-SO 

ATTN 

AUL-LSE 

US  Army  Nuclear  A  Chemical  Agency 

Deputy  Chief  of  Staff 

ATTN:  Library 

Rsch,  Dev  A  Acq 

3  cy  ATTN 

AFRDS,  Space  Sys  A  C3 

US  Army  Research  Ofc 

ATTN:  R.  Mace 

Rome  Air  Dev  Ctr 

ATTN 

0C0,  J.  Simons 

DEPARTMENT  OF  THE  NAVY 

Space  Div 

Naval  Electronics  Systems  Cmd 

ATTN 

WE 

ATTN:  PNE  117-20 

ATTN 

YGD 

ATTN 

YN,  P.  Sivgals 

Naval  Intelligence  Support  Ctr 

ATTN:  Doc  Con 
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DEPARTMENT  OF  ENERGY 

Department  of  Energy 

Office  of  Military  Application 
ATTN:  OMA,  DP-22 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 

Department  of  Commerce 

National  Bureau  of  Standards 
ATTN:  A.  Phelps 

Department  of  Commerce 

National  Bureau  of  Standards 

ATTN:  Sec  Ofc  for  S.  Abramowitz 
ATTN:  Sec  Ofc  for  R.  Levine 
ATTN:  Sec  Ofc  for  M,  Krauss 
ATTN:  Sec  Ofc  for  J.  Devoe 

Department  of  Commerce 

National  Oceanic  &  Atmospheric  Admin 
3  cy  ATTN:  F.  Fehsenfeld 
3  cy  ATTN:  E.  Ferguson 


DEPARTMENT  OF  ENERGY  CONTRACTORS  (Continued) 

Los  Alamos  National  Lab 

ATTN:  MS670,  J.  Hopkins 
ATTN:  Librarian 
ATTN:  G.  Smith 
ATTN:  MS362,  Library 
ATTN:  R.  Jeffries 

Sandia  National  Labs 
ATTN:  T.  Cook 

Sandia  National  Labs 

ATTN:  Tech  Library,  3141 
ATTN:  M.  Kramm 
ATTN:  ORG  1250,  W.  Brown 
ATTN:  L.  Anderson 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aero-Chem  Research  Labs,  Inc 
ATTN:  A.  Fontijn 

Aerodyne  Research,  Inc 
ATTN:  C.  Kolb 
ATTN:  M.  Camac 


Institute  for  Telecommunications  Sciences 


Aerospace  Corp 


ATTN:  W.  Utlaut 

ATTN 

N.  Cohen 

ATTN:  G.  Falcon 

ATTN 

V.  Josephson 

ATTN 

I.  Garfunkel 

NASA 

ATTN 

J.  Reinheimer 

ATTN:  Code  6801,  A.  Tempkin 

ATTN 

J.  Straus 

ATTN:  Tech  Library 

ATTN 

T.  Widhoph 

ATTN:  Code  900,  J.  Siry 

ATTN 

H.  Mayer 

3  cy  ATTN:  A.  Aiken 

AYCO  Everett  Research  Lab,  Inc 

NASA 

ATTN 

mo 

ATTN:  N.  Stone 

ATTN 

C.  Von  Rosenberg,  Jr 

ATTN:  W.  Roberts 

ATTN 

Tech  Library 

ATTN:  J.  Watts 

Batelle  Memorial  Institute 

NASA 

ATTN 

H.  Lamuth 

ATTN:  MS-168,  C.  Schexnayder 

ATTN 

Stoiac 

ATTN 

R.  Thatcher 

NASA 

ATTN:  N-245-3,  R.  Whitten 

Berkeley  Rsch  Assoc,  Inc 

ATTN 

J.  Workman 

NASA 

ATTN:  J.  Gray 

Boston  College 

ATTN 

E.  Heglom 

NASA  Headquarters 

ATTN 

W.  Greider 

ATTN:  I.  Schardt.  Code  EE 

The  Trustees  of  Boston  College 

OTHER  AGENCY 

2  cy  ATTN 

Chairman  Dept  of  Chemisti 

2  cy  ATTN 

Chairman  Dept  of  Physics 

Government  Publications  Library-M 

ATTN:  J.  Winkler 

University  of  California  at  Riverside 

ATTN:  J.  Pitts,  Jr 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Cal  Ifornla 

Institute  of  Technology 

University  of  California 

ATTN:  J.  Ajello 

Lawrence  Livermore  National  Lab 

ATTN:  L-71,  J.  Chang 

Cal  span  Corp 

ATTN:  L-262,  0.  Wuebbles 

ATTN:  M.  Dunn 

ATTN:  L-10,  H.  Kruger 

ATTN 

J.  Grace 

ATTN:  L-48,  E.  Woodward 

ATTN 

C.  Treanor 

ATTN:  L-325,  G.  Haugan 

ATTN 

W.  Wurster 
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DEPARTMENT  OF  OEFENSE  CONTRACTORS  (Continued) 

University  of  Colorado 

ATTN:  C.  Lineberger-JILA 
ATTN:  G.  Lawrence-LASP 

Columbia  University 

ATTN:  Security  Officer  for  H.  Foley 

Computer  Sciences  Corp 

ATTN:  F.  Eisenbarth 

Concord  Sciences 

ATTN:  E.  Sutton 

Cornell  University 

ATTN:  M.  Kelly 

University  of  Denver 

ATTN:  Sec  Officer  for  D.  Murcray 

University  of  Denver 

ATTN:  B.  Van  Zyl 

Environmental  Rsch  Inst  of  Michigan 
ATTN:  IRIA,  library 

EOS  Technologies,  Inc 

ATTN:  B.  Gabbard 

General  Electric  Co 

ATTN:  R.  Edsall 
ATTN:  P.  Zavitsanos 

General  Research  Corp 
ATTN:  J.  Ise,  Jr 

Geo  Ctrs,  Inc 

ATTN:  E.  Marram 

Institute  for  Defense  Analyses 
ATTN:  H.  Wo If hard 
ATTN:  E.  Bauer 

Kaman  Sciences  Corp 

ATTN:  J.  Jordano 

Kaman  Sciences  Corp 

ATTN:  T.  Stephens 


Kaman  Tempo 
ATTN: 
ATTN: 
ATTN: 
5  cy  ATTN: 


B.  Gambill 

J.  Devore 

K.  Schwartz 
DAS  I AC 


Kaman  Tempo 

ATTN:  DASIAC 

Lockheed  Missilies  &  Space  Co,  Inc 
ATTN:  B.  McCormac 
ATTN:  M.  Walt 
ATTN:  J.  Cladis 
ATTN:  J.  Reagan 
ATTN:  J.  Perez 
ATTN:  R.  Sears 
4  cy  ATTN:  J.  Kumer 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

University  of  Lowell 
ATTN:  G.  Best 

Mission  Research  Corp 


ATTN 
ATTN 
ATTN 
ATTN 
5  cy  ATTN 


M.  Scheibe 

O,  Archer 

F.  Guigliano 

P.  Fischer 
Tech  Library 


Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  D.  Davis 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

PhotoMetrics,  Inc 

ATTN:  I.  Kofsky 

Physical  Dynamics,  Inc 
ATTN:  A.  Thompson 

Physical  Science  Lab 

ATTN:  W.  Berning 

Physical  Sciences,  Inc 
ATTN:  R.  Taylor 
ATTN:  K.  Wray 
ATTN:  G.  Calendonia 

Physics  International  Co 
ATTN:  Tech  Library 

University  of  the  Commonwealth,  Pittsburgh 
ATTN:  F.  Kaufman 
ATTN:  M.  Biondi 
ATTN:  W.  Fite 

Princeton  University 
ATTN:  Librarian 

R  !  D  Associates 

ATTN:  F.  Gilmore 
ATTN:  R.  Turco 
ATTN:  R.  Lindgren 
ATTN:  H.  Ory 
ATTN:  P.  Haas 

R  &  D  Associates 

ATTN:  B.  Yoon 
ATTN:  J.  Rosengren 

Rand  Corp 

ATTN:  P.  Davis 
ATTN:  C.  Crain 

Rand  Corp 

ATTN:  B.  Bennett 

Science  Applications,  Inc 
ATTN:  D.  Sachs 
ATTN:  D.  Hamlin 

Space  Data  Corp 

ATTN:  S.  Fisher 


1 


1 


DEPARTMENT  OF  OEFENSE  CONTRACTORS  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Spectral  Sciences,  Inc 
ATTN:  F.  Beln 


Technology  international  Corp 
ATTN:  W.  Boqulst 


SRI  International 


Utah  State  University 


ATTN 

D.  McDaniels 

Attention  Sec  Con  Ofc  for 

ATTN 

W.  Chesnut 

ATTN 

C. 

Wyatt 

ATTN 

M.  Baron 

ATTN 

0. 

Baker 

ATTN 

R.  Leada brand 

ATTN 

A. 

Steed 

ATTN 

J.  Casper 

ATTN 

K. 

Baker,  Dir 

SRI  International 

IfisiDyne,  Inc 

ATTN:  C.  Hulbert 

ATTN 

W. 

Reidy 

ATTN 

C. 

Humphrey 

•Stewart  Radiance  Lab 

ATTN 

T. 

Degges 

ATTN 

R.  Huppi 

ATTN 

H. 

Smith 

ATTN 

J. 

Carpenter 

Teledyne  Brown  Engrg 

ATTN 

MS-12,  Tech  Library 

Wayne  State  University 

ATTN 

F.  Leopard 

ATTN 

R. 

Kunmler 

ATTN 

N.  Passino 

t  Space  Sci 


Lockheed  Missiles  &  Space  Co,  Inc 
4  cy  ATTN;  T.  James 


Wayne  State  University 
ATTN:  W.  Kauppila 


